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Semi-Annual  Progress  Report 

THREE  DIMENSIONAL  TRANSONIC  Fl.OWS  (N  COMPRESSORS  AND  CHANNELS 

The  University  of  Michigan,  Ann  Arbor,  Micliigan 
Subcontract  No.  8960-10 


Professor  T.  C.  Adamson,  Jr.,  Principal  Investigator 
Professor  M.  Sichel,  Principal  Investigator 

Introd uctlo n 

The  work  described  here  is  a study  of  three  dimensional  transonic  chan- 
nel flows  using  asymptotic  methods.  Consideration  of  such  flows,  which  a- 
rise  in  transonic  cascades,  is  Important  because  of  tlie  possibility  of  mixed 
supersonic-subsonic  flows  where  two  dimensional  solutions  are  completely  In- 
adequate. The  analysis  has  focused  on  those  regimes  of  flow  where  asymptotic 
methods  can  be  used  to  derive  relatively  simple  analytical  solutions  to  the 
flow  fields.  A model  flow  consisting  of  a shear  flow  passing  through  a three 
dimensional  channel  with  a flow  constriction  has  been  considered.  This  model 
problem  contains  the  main  features  of  flow  in  a cascade,  but  in  a very  simple 
geometry.  Asymptotic  methods  have  been  used  to  identify  those  regions  where 
solutions  derived  from  linearized  equations  are  not  uniformly  valid,  and  where, 
therefore,  it  becomes  necessary  to  use  nonlinear  equations.  The  basic  philo- 
sophy, has  been  to  develop  simple  approximate,  but  analytical,  solutions  which 
display  the  effects  of  various  parameters  on  the  nature  of  the  flow,  rather 
than  to  attempt  a detailed  numerical  analysis. 
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Discussion 


The  model  problem  chosen  for  study  may  be  interpreted  as  the  flow  through 
a linear  cascade  with  the  blades  aligned  parallel  to  the  incoming  flow.  The 
symmetry  boundaries  upstream  and  downstream  of  the  blades  are  replaced  by  walls 
so  that  the  flow  considered  is  that  through  a rectangular  channel  with  flow 
constrictions  corresponding  to  half  blades  on  opposite  walls.  The  radial  vari- 
ation of  the  rotor  tangential  velocity  component  is  represented  by  a linear 
gradient  of  the  velocity  entering  the  channel.  Because  of  this  gradient,  a 
portion  of  the  incoming  flow  may  be  subsonic  and  a portion  supersonic,  and  this 
is  the  type  of  flow  which  has  been  considered. 

The  nature  of  this  model  flow  will  depend  on  the  location  of  the  sonic 
surface  of  the  incoming  shear  flow  with  respect  to  the  channel,  the  thickness 
ratio  of  the  flow  constrictions  used  to  represent  the  turbine  blades,  and  the 
velocity  gradient  of  the  incoming  shear  flow.  The  essential  features  of  such 
flows  in  the  absence  of  shock  waves  are  described  in  references  1 and  2.  Key 
results  are  the  verification  tliat  linearized  equations  can  be  used  to  des- 
cribe most  of  the  flow  and  the  establishment  of  conditions  under  which  chok- 
ing occurs.  In  addition,  the  inner  regions  at  the  blade  leading  edges,  where 
weak  shocks  may  occur,  and  at  the  plane  of  minimum  crosssectional  area,  where 
singularities  in  the  solution  may  occur,  have  been  investigated.  These  re- 
sults are  described  in  references  1 and  2. 

Shock  waves  may  appear  when  a flow  with  a subsonic  average  Mach  number 
at  the  inlet  accelerates  until  an  average  supersonic  Macii  number  is  attained 
dowiiMlream  of  the  si-cllon  of  minimum  area.  Tlie  detailed  study  of  such  shock 
wavi's  has  been  continued  during  tlie  period  covered  by  tills  progress  report. 

If  the  flow  accelerates  to  a slightly  supersonic  average  Mach  numher,  the 
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shock  wave  which  may  appear  will  not  fill  the  channel  because  of  the  gradient 
in  the  incoming  flow.  This  feature  makes  the  analysis  of  such  shocks  more 
complicated  than  in  the  simpler  case  of  one  dimensional  channel  flow. 

As  established  earlier,  most  of  the  flow  field  can  be  described  by  linear 
equations;  but,  very  close  to  the  shock  wave,  the  nonlinear  transonic  small 
disturbance  equations  are  required.  In  order  to  understand  the  main  general 
features  of  the  flow,  the  method  of  Integral  relations  is  being  used  to  de- 
velop simple  approximate  solutions.  Preliminary  results,  which  have  been 
obtained  so  far,  show  that  if  the  sonic  line  reaches  a position  sucli  that 
90%  of  the  channel  flow  is  supersonic  just  ahead  of  the  sliock,  then  the  siiock 
extends  about  75%  of  the  distance  across  the  channel.  On  the  other  hand, 
tiie  shock  extends  only  about  32%  of  the  distance  across  the  channel  when  70% 
of  the  mainstream  flow  is  supersonic.  The  detailed  analysis  leading  to  these 
results  have  revealed  the  complexity  of  such  flows,  particularly  with  respect 
to  the  region  of  upstream  Influence  of  the  flow  downstream  of  the  shock. 

As  already  mentioned,  there  is  a region  near  the  leading  edge  of  the 
blades  where  the  linear  equations  used  In  most  of  the  flow  arc  invalid.  Ttie 
flow  in  this  region  has  been  studied  in  detail  during  tlie  current  contract 
period  and  it  has  been  possible  to  show  that  local  solutions  at  the  leading 
edge  match  asymptotically  with  the  first  order  solutions,  used  in  the  earlier 
analysis,  which  are  valid  away  from  the  edge. 

References 

(1)  Adamson,  Jr.,  T.  C.,  "Three  Dimensional  Transonic  Shear  Flow  in  a Chan- 
nel," Transonic  Flow  Problems  in  Turbomachinery,  eds.  T.  C.  Adamson, 

Jr.  and  M.  F.  Platzer,  Hemisphere  Publishing  Corporation,  1977,  p.  70. 

(2)  Slchel,  M. , and  Adamson,  Jr.,  T.  C. , "Three  Dimensional  Transonic  Flow 
in  Channels,"  XIII  Symposium  on  Advanced  Problems  and  Methods  in  Fluid 
Mechanics,  Olsztyn-Kortowo,  September  5-10,  1977. 
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AXIAL  FLOW  FAN  S'I'ACL  UNSTKADY  I’F RFORMANCF 


Applied  Research  Laboratory 
The  Pennsylvania  State  University 
P.  0.  Box  30,  State  College,  Pennsylvania  16801 

Subcontract  No.  8960-A 

Edgar  P.  Bruce,  Principal  Investigator 


Introduction 

The  objective  of  this  research  is  to  analyze  the  time-dependent 
interaction  between  the  components  of  an  Isolated  axial  flow  fan  stage 
and  a spatially  fixed,  circumferentially  varying  flow  field.  The  major 
variables  are  reduced  frequency;  rotor  blade  space-to-chord  ratio,  stagger 
angle,  mean  angle  of  attack,  and  design  loading  level;  and  rotor-stator 
axial  spacing. 

The  experiments  are  being  conducted  in  the  ARL  Axial  Flow  Research 
Fan.  This  facility  has  a hub  radius  of  12.06  cm  (A. 75  inches),  a hub-to-tlp 
radius  ratio  of  0.AA2,  and  operates  in  the  subsonic  incompressible  flow  regime. 
The  rotor  and  stator  blades  have  a 10  percent  thick  Cl  profile  with  a chord 
of  15. 2A  cm  (6.00  inches)  and  an  aspect  ratio  of  unity. 

Instrumentation  available  at  present  or  under  development  consists  of; 

(1)  a strain  gaged  sensor  mounted  within  one  rotor  blade  which  detects  the 
time-dependent  normal  force  and  pitching  moment  developed  on  a mid-span 
blade  segment,  (2)  hot-film  sensors  mounted  on  the  suction  surface  of 
rotor  and  stator  blades  which  detect  the  nature  of  the  boundary  layer, 
i.e.,  whether  the  Instantaneous  boundary  layer  flow  is  laminar,  turbulent 
or  separated;  (3)  dynamic  total  head  probes;  (A)  two-element  hot-film 
probes;  and  (5)  conventional  three-dimensional  directional  probes.  A 
system  is  being  developed  which  will  permit  on-line  analysis  of  all  time- 
dependent  signals  by  a digitizing,  phase-lock  averaging  process. 

The  unsteady  normal  force  and  pitching  moment  results  obtained  in  the 
initial  phase  of  this  program  at  reduced  frequencies  from  0.2  to  2.1  have 
been  documented  in  a Project  SQUID  report  (Reference  1) . Since  completing 
the  initial  phase,  our  efforts  have  been  directed  toward  extending  the 
reduced  frequency  range  of  the  uncambered  rotor  experiments  (Reference  1) 
from  2.1  to  5.0,  and  toward  a detailed  examination  of  the  effects  of  inflow 
distortion  on  the  performance  of  a stage  designed  with  a free-vortex  loading 
distribution.  A related  theoretical  effort  has  as  its  goal  an  extension  of 
the  unsteady  lift  cascade  model  developed  by  Henderson  (Reference  2)  to 
include  the  unsteady  pitching  moment. 
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Dlscusslon 

During  this  reporting  period,  we  have  made  significant  progress  in  both 
the  experimental  and  theoretical  areas.  With  respect  to  the  experimental  tasks, 
we  have: 

1)  completed  through  final  data  reduction  a thorough  recalibration 
of  all  the  disturbance  generating  screens  used  to  produce  multi- 
cycle rotor  inlet  axial  velocity  profiles, 

2)  measured,  using  5-hole  probes  located  immediately  upstream  of 
and  approximately  one-half  chord  length  downstream  of  the  free- 
vortex  loading  distribution  rotor,  the  circumferential  variation 
of  the  time  mean  inlet  and  exit  flow  properties  at  mean  incidence 
angles  of  zero  and  seven  degrees  with  different  disturbance 
generating  screens  installed  to  give  a reduced  frequency  range 
from  0 to  5.0,  and 

3)  completed  static  and  dynamic  calibration  of  the  sensors  required 
for  the  unsteady  force  and  moment  testing  of;  a)  the  free-vortex 
loading  distribution  rotor  at  reduced  frequencies  from  0 to  5.0, 
and  b)  the  uncambered  rotor  at  reduced  frequencies  from  2.1  to 

5.0. 

Data  from  the  testing  outlined  under  2)  above  are  being  reduced  at  present. 
Preliminary  results  Indicate  that  a reduced  frequency  range  exists  for  which 
an  incoming  sinusoidal  distortion  in  axial  velocity  component/total  pressure 
is  attenuated  to  approximately  one-tenth  of  its  upstream  level.  On  either  side 
of  this  reduced  frequency  range,  the  corresponding  flow  properties  are 
attenuated  to  approximately  forty  percent  of  their  upstream  levels.  The 
testing  outlined  under  3a)  above  will  be  initiated  in  March  1978  and  will  be 
immediately  followed  by  the  testing  outlined  under  3b). 

An  independent  derivation  of  the  expression  developed  by  Henderson 
(Reference  2)  for  the  unsteady  lift  generated  on  the  blades  of  a cascade 
due  to  operation  in  a distorted  inflow  has  been  completed  as  a prelude  to 
extension  of  this  model  to  include  an  expression  for  the  unsteady  pitching 
moment.  The  results  of  the  present  derivation,  in  general,  agree  with 
Henderson's  results;  however,  errors  have  been  found  and  corrected  in  the 
following  areas: 

1)  the  combination  of  values  of  the  cascade  parameters  space-to-chord 
ratio  (S/C)  and  stagger  angle  (C)  for  which  the  results  are  valid. 
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2)  the  derivation  of  one  of  the  two  terms  which  define  the  effect 
on  the  unsteady  lift  on  the  reference  blade  (n  = 0)  due  to  all 
the  neighboring  blades  (-  “>  < n < - 1)  and  1 < n < <») , and 

3)  the  type  of  camber  line  for  which  the  solution  is  valid. 

The  effect  of  the  change  und^r  item  1)  above  is  to  include  certain  values  of 
(S/C)  that  were  previously  excluded  for  values  of  sin  C < 0.25  and  to  exclude 
certain  values  of  (S/C)  that  were  previously  Included  for  values  of  sin  > 0.25. 
A preliminary  evaluation  of  the  effect  on  predicted  unsteady  lift  due  to  the 
change  noted  under  item  2)  above  indicates  that  the  effect  is  small.  With 
regard  to  item  3)  above,  the  camber  line  shape  employed  in  Reference  2 is  a 
symmetrical  parabolic  arc  — not  a circular  arc. 

The  unsteady  moment  expression  has  been  derived  to  the  point  where  the 
expression  for  the  symmetrical  parabolic  arc  camber  line  shape  and  the  flow 
tangency  boundary  condition  are  being  included.  After  simplifying  the  resulting 
expression  as  much  as  is  possible,  it  will  be  programmed  for  solution  on  the 
IBM  360  computer. 


^^f^ren^^s 

1.  Bruce,  E.  P.  and  Henderson,  R.  E.,  "Axial  Flow  Rotor  Unsteady 
Response  to  Circumferential  Inflow  Distortions,"  Project  SQUID 
Technical  Report  PSU-13-P,  September  1975. 

2.  Henderson,  R.  E.,  "The  Unsteady  Response  of  an  Axial  Flow  Turbo- 
machine to  an  Upstream  Disturbance, ", PI). D.  Dissertation,  Engineering 
Department,  University  of  Cambridge,  1972. 
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INVESTIGATION  OF  THE  EFFECTS  OF  HIGH  AERODYNAMIC 
LOADING  ON  A CASCADE  OF  OSCILLATING  AIRFOILS 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06] OR 

I Subcontract  896O-19 

Franklin  0.  Cai'ta,  Principal  Inventifp-itor 
Arthur  0.  St.  Hilaire,  Principal  Investigator 

Introduction 

The  basic  objective  of  this  research  program  is  to  study 
the  phenomenon  of  dynamic  stall  on  a cascade  of  oscillating 
airfoils.  Measurements  are  being  made  of  the  unsteady  chordwise 
pressure  distribution,  and  efforts  are  being  made  to  detect 
the  occurrence  of  boundary  layer  transition  and  separation 
on  the  surface  of  an  oscillating  cascaded  airfoil  operating 
near  the  static  stall  condition. 

Program  Review 

-O 

Tents  completed  at  a mesin  incidence  angle  of  b continue  to 
show  the  strong  influence  of  interblade  phase  angle  on  blade 
aerodynamic  response.  This  is  based  on  oscilloscope  observations 
made  during  the  test;  however,  it  is  evident  that  the  stall 
breakdown  of  the  pressure  signal  near  the  leading  edge  is  altered 
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radically  as  interblade  phase  anple  is  chanp:ed  from  positive  to 
negative  values. 

In  addition  to  blade  unsteady  pressures , measurements  are  also 
being  made  of  the  unsteady  sidewall  pressures  in  the  plane  of  the 
cascade.  Strong  oscillatory  signals  have  been  observed  in  both 
leading  and  trailing  edge  regions.  These  data  are  needed  to 
assist  in  the  formulation  of  an  analytical  model  for  unsteacly 
flows  past  heavily  loaded  cascades,  which  is  the  basis  of  our 
planned  future  work. 

Upon  completion  of  our  current  test  at  10  incidence,  the 
system  will  be  prepared  for  our  phase  11  work  vdiich  consists 
of  selecting  a specific  mean  incidence  angle,  testing  at 
several  gap/chord  ratios,  and  reducing  all  of  the  data. 
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INVESTKiATlON  OF  ADVERSE  PRESSURE  GRADIENT  CORNER  FLOWS 


University  of  Washington,  Seattle,  Washington 
Subcontract  No.  8960-27 


Professor  F.  B.  Gessner,  Principal  Investigator 
Mr.  S.  Ono,  Research  Assistant 


Introduct ion 

This  program  is  a comprehensive  experimental  study  of  turbulent  flow  along 
a St  Teamwise  corner  in  the  presence  of  an  adverse  jiressure  gradient.  Under  those 
conditions  it  is  known  that  secondary  flows  appear  in  the  vicinity  of  the  corner 
and  arc  amplified  as  the  flow  decelerates  [1].  For  moderate  adverse  pressure 
gradient  (non-separating)  flow  conditions,  secondary  flow  velocities  arc  approxi- 
mately 10  per  cent  of  the  primary  flow  component,  which  is  an  order  of  magnitude 
greater  than  that  which  has  been  observed  in  favorable  and  zero  pressure  gradient 
corner  flows.  For  conditions  near  incipent  separation,  it  is  probable  that  the 
10  per  cent  figure  noted  above  will  be  exceeded,  especially  in  the  vicinity  of 
a corner  where  the  flow  first  separates  locally  from  the  wall. 

In  order  to  make  measurements  in  an  environment  where  the  mean  flow  direction 
varies  from  point  to  point  in  the  flow,  one  must  decide  whether  to  align  probes 
with  the  local  mean  flow  direction  or  account  for  cross-flow  effects  using  probes 
aligned  with  the  streamwise  flow  component.  The  latter  method  of  measurement  was 
chosen  in  the  present  study,  primarily  because  of  probe  alignment  difficulties 
associated  with  the  former  method.  In  order  to  make  accurate  measurements  in  a 
skewed  mean  flow,  hot-wire  probe  configurations  must  be  selected  which  minimize 
distortion  of  flow  in  the  vicinity  of  the  sensing  element  induced  by  flow  over 
the  supports.  The  pitch  and  yaw  angles  of  the  wire  relative  to  the  local  mean 
flow  direction  must  also  be  such  that  cross-flow  induced  tangential  cooling 
effects  can  be  readily  modelled.  In  the  present  study  a single  "normal"  wire 
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probe  (a=  0°)  and  a rotatable  inclined  wire  probe  (a=  45°)  were  selected  in 
lieu  of  an  X array  or  a multi-wire  configuration  to  minimize  drift  problems  and 
resolution  inaccuracies.* 


Discussion 


For  the  past  six  months  we  have  been  engaged  in  an  extensive  experimental 
program  designed  to  provide  information  on  response  behavior  of  hot-wires  in  a 
skewed  mean  flow.  The  extent  of  this  work  exceeded  our  original  expectations, 
but  has  proven  to  be  we  1 1 -worthwhile  from  the  standpoint  of  providing  vital 
information  needed  to  make  accurate  mean-flow  and  turbulence  measurements  in 
this  environment.  In  general,  the  mean  bridge  voltage  output  from  a hot-wire 
anemometer,  lij^,  may  be  expressed  ns 


or,  alternatively,  as 


= r;^  + BU“'^  + CU 

b o c e 


(2) 


where  is  the  effective  cooling  velocity  across  the  wire,  n is  an  exponent 

which  varies  typically  from  0.4  and  0.5,  B and  C arc  coefficients  which  arc 

constant  for  constant  temperature  operation  of  the  wire,  and  is  the  inter- 

cept  of  a Icast-squares-fi t of  the  data  when  plotted  either  as  vs  or  as 

vs  u'’’^  . The  effective  cooling  velocity  can  be  modelled  in  terms  of  the 
be  b ^ 

local  mean  velocity  component  in  a given  plane,  U^,  the  angle  between  and 
the  normal  to  a wire  in  that  plane,  iji  , and  a tangential  cooling  factor,  k,  as 
fol lows: 

= U^(cos^  ij)  + k^  sin^  ip)'^  (.^) 

For  turbulence  measurements  with  inclined  wires  it  is  imperative  that  the 
factor  k be  determined  accurately  for  each  local  operating  condition,  inasmuch 
as  turbulent  shear  stress  components  arc  sensitive  to  the  ratio  (1+  k^)/(l-k^) 

*Thc  angle  a is  defined  as  the  angle  between  the  streamwisc  component  of  the 
local  mean  velocity  in  a given  plane  and  the  normal  to  a wire  in  that  plane. 
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and  transverse  normal  stress  components  are  sensitive  to  ( l+k^) ^/ ( 1 -k^ ) ^ when 
a=  45°  [2].  In  applications  where  an  inclined  wire  is  oriented  to  i|;  = u = 45° 

(no  cross  flow),  it  is  permissablc  to  assume  that  k=  k(f/d)  only  over  a restricted 
range  of  values  for  U^,  where  £/d  is  the  wire  Icngth-to-d iameter  ratio  [51.  When 
the  mean  flow  is  skewed,  however,  k is  also  deiiendent  on  i|),  so  tliat  k=-  kCf/d.ili)  I II. 
These  variations  of  k with  tp  arc  not  small,  and  can  leatl  to  a 2:1  change  in  k values 
as  i|/  varies  from  15  to  75  deg-^ees.  In  References  .5  and  4 only  relatively  narrow 
ranges  of  velocity  were  considered,  and  a constant  value  of  n(0.45)  was  assumed 
to  be  adequate  for  correlating  the  data  in  each  velocity  range. 

In  our  work  we  have  examined  the  scope  of  applicability  of  these  correlations 
over  a relatively  wide  range  of  flow  conditions  (1  _<  £ 30  mps,  15  < p £ 75°). 

This  work  entailed  modification  of  our  jet  flow  calibration  facility  by  the  addition 
of  a viewing  microscope  and  turntable.  With  these  modifications  wire  inclination 
angles  could  be  set  with  precision  to  within  ±0.25  degrees.  Our  turbulent  pipe 
flow  calibration  facility  was  also  modified  by  adding  a traversing  and  rotating 
mechanism  at  the  pipe  discharge.  This  modification  enabled  us  to  simulate  a 
skewed  mean  flow  across  a wire  by  inclining  the  probe  in  fixed  degree  increments 
from  the  axial  flow  direction. 

With  5 p diameter  tungsten  wire  probes  located  in  the  potential  core  of  a 
free  jet  (170  ^/d  £ 233),  the  moan  bridge  voltage  was  measured  for  various 

combinations  of  and  ip.  A least-squares-fit  was  applied  to  the  data  and  k 
was  evaluated  by  means  of  either  equation  (1)  or  equation  (2)  in  combination 
with  equation  (3) . The  results  indicate  that  when  equations  (1)  and  (3)  are 
utilized,  k(i|))  is  a uniquely  specified  function,  provided  that  the  variation 
of  n with  is  prescribed  properly.  The  exponent  n varies  typically  from 
n=0.48  at  U^=  3 mps  to  n~0.40  when  U^=  30  mps.  These  results  have  important 
implications  because  they  demonstrate  that  indiscriminate  use  of  the  k(£/d) 
correlation  presented  by  Champagne  et  al.[3]  and  the  k(£/d,ij;)  correlation  devel- 
oped by  Friehe  and  Schwarz  [4]  will  lead  to  significant  errors  in  estimated  values 
for  k unless  n is  allowed  to  vary  with  U^,  even  when  no  cross  flow  exists. 

The  adequacy  of  the  correlations  for  k(£/d,ij;)  and  n(U^)  which  have  been 
developed  in  the  present  study  have  been  tested  by  means  of  measurements  with 
a rotatable  inclined  wire  probe  in  the  potential  core  of  a free  jet.  The  results 
indicate  that  pitch  and  yaw  angles  as  large  as  30  degrees  can  be  measured  to 
within  an  accuracy  of  ±1  degree  when  the  prescribed  correlations  are  used.  If 
a mean  (constant)  value  for  k is  specified,  pitch  and  yaw  angles  are  overestimated 
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by  as  much  as  5 degrees.  In  order  to  analyze  the  response  sensitivity  of 
hot-wires  to  turbulent  fluctuations  in  a skewed  mean  flow,  response  equations 
(second-order  accurate)  have  been  developed  which  account  for  tangential  cooling 
effects  induced  by  a transverse  mean  flow.  The  range  of  validity  of  these 
equations  was  examined  by  making  Reynolds  stress  measurements  in  fully  developed 
turbulent  pipe  flow  with  "normal"  and  inclined  wires  oriented  at  yaw  angles 
of  10,  20,  and  30  degrees  relative  to  the  axial  mean  flow  direction.  The  measure- 
ments were  made  over  a radial  interval  0.1  _<  r/R  < 1.0  at  bulk  flow  Reynolds 
numbers  of  20,000  and  150,000  for  which  n=^  0.48  and  n=^  0.40,  respectively.  The 
results  show  that  first-order  theory  begins  to  break  down  when  yaw  angles  exceed 
10  degrees,  and  that  second-order  theory  is  accurate  for  yaw  angles  as  high  as 
30  degrees,  provided  that  kff/d,  ij^)  and  nCU^)  are  specified  in  accordance  with 
the  correlations  whicli  have  been  developed.  This  phase  of  the  work  will  be 
written  u])  shortly  in  the  form  of  a Project  SQUID  Technical  Report. 

The  above  study  has  delayed  somewhat  the  initiation  of  measurements  in 
oiir  rectangular  diffiuer  facility  for  the  purpose  of  investigating  adverse 
Iircssure  gradient  effects  on  flow  in  the  corner  region.  Wc  feel  justified  in 
pursuing  our  response  sensitivity  study,  however,  because:  (1)  the  work  is  of 
fundamental  importance,  and  (2)  the  information  gleaned  from  this  phase  of  the 
work  will  assist  us  in  analyzing  data  obtained  in  our  corner  flow  work. 
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Introduction 


The  general  goal  of  this  program  is  to  study  the  transitory  stall 
flow  regime  in  two-dimensional  diffusers.  Maximum  value  of  pressure 
recovery  at  fixed  non-dimensional  length,  an  important  design  optimum 
[1],  generally  occurs  when  the  turbulent  boundary  layers  are  starting  to 
separate  or  stall.  The  flow  is  rather  unsteady  and  significant  amounts 
of  transient  back  flow  already  are  seen  in  the  diffuser  at  peak  pressure 
recovery.  These  flow  conditions  are  associated  with  the  onset  and 
development  of  the  transitory  stall  flow  regime  [2]. 

Ghose  and  Kline  [3]  have  developed  a new,  steady  flow  boundary 
layer  prediction  method  which  is  solved  simultaneously  (not  iteratively) 
with  the  inviscid  core  flow.  This  method  gives  surprisingly  good  agree- 
ment with  data  on  pressure  recovery  up  to,  and  slightly  beyond  the 
condition  of  peak  recovery.  The  existing  wall  pressure  data  in  this 
region  are  not  of  sufficient  accuracy  to  properly  check  the  method, 
however. 

The  primary  objectives  of  our  program  are  (i)  to  provide  new  mean 
and  fluctuation  velocity  and  pressure  data  in  diffusers  operating  close 
to  peak  pressure  recovery  in  order  to  complement,  check  and  provide  a 
data  base  of  sufficient  accuracy  to  allow  for  possible  improvement  of 
the  prediction  method  of  Ghose  and  Kline  [3],  and  (ii)  to  study  the 
magnitude  of  the  velocity  and  pressure  fluctuations  in  the  transitory 
stall  regime  in  order  to  provide  a useful  extension  of  the  work  of  Smith 
and  Kline  [2]  and  Layne  and  Smith  [4], 
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Discussion 


Work  is  proceeding  in  several  areas,  (i)  the  collection  and 
correlation  of  wall  static  pressure  distributions  and  inlet  velocity 
profiles,  and  (ii)  the  development  of  a new  wall  flow  direction  instru- 
ment. 


The  diffuser  tunnel  is  now  complete  except  for  the  heat  exchanger 
and  associated  control  components.  It  has  been  checked  out  and  is  now 
being  used  to  collect  the  mean  flow  data. 

Diffuser  Experiments  are  now  in  progress  with  the  objective  of  obtain- 
ing accurate,  detailed  static  pressure  distributions  in  a long 
(N/W^  = 15)  , straight-walled  diffuser  for  a set  of  opening  angles,  20  . 

Fig.  1 in  the  progress  report  of  March  15,  1977,  shows  the  general 
configuration  for  this  first  series  of  tests.  The  diffuser  inlet  width 
is  W]  = 3.00  inches  and  the  distance  between  parallel  walls  is 
b = 12.00  inches  so  the  inlet  aspect  ratio  is  4:1  . 

Fig.  1,  below,  shows  some  details  of  the  inlet  configuration.  The 
13.4  inch  long  inlet  duct,  with  its  boundary  layer  trips,  assures  the 
development  of  identical  turbulent  boundary  layers  on  all  four  diffuser 
walls  at  the  inlet  station  located  a distance  S/W]  = -1.15  with  respect 
to  the  diffuser  throat. 

For  the  current  series  of  tests,  the  measured  inlet  conditions  are: 

Air  Pressure:  1 atmosphere 
Air  Temperature:  84°F  ^ 

Kinematic  Viscosity:  1.7  x 10"^  ft'^/s 
Free  Stream  Speed:  Ug]  = 151  ft/s 
Reynolds  Number:  W-|Ugi/v  = 2.2  x 10^ 

Boundary  Layer  Parameters  (Preliminary): 

6 (99X)  = 0.33  inches 
6*  = 0.041  inches 
0 = 0.032  inches 

H = 1.3 

The  overall  static  pressure  recovery  coefficients,  CpyQ.j.  = (Pi-P2)/*5pU  , , 
are  given  below.  The  inlet  state  (1)  is  defined  at  the  ® 

inlet  station  denoted  on  Fig.  1.  The  exit  station  (2)  is  located  in 
the  short,  constant  area  tailpipe,  3 inches  downstream  of  the  diffuser. 
Experimental  uncertainties  for  the  overall  recovery  coefficients  are  yet 
to  be  determined,  but  we  hope  to  achieve  uncertainties  as  low  as  + 1 to 
2t  (20:1  odds). 
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2()  (degrees) 

^PfOT 

Remarks 

0 

— 

Test  to  establish  inlet  pressure  drop 

4 

Unstalled  case 

8 

.714 

Early  transitory  stall 

10 

— 

Close  to  peak  Cp^Qj 

12 

.655 

Transitory  stall 

16 

.610 

Transitory  stall  with  large  fluctuations 

20 

.538 

Transitory  stall  with  large  fluctuations 

24 

— 

Transitory  stall  with  large  fluctuations 

The  detailed,  local  wall  static  pressure  distribution  for  the  results 
obtained  to  date  are  shown  in  Fig.  2.  Work  is  in  progress  to  complete 
the  results  at  20  = 0,  4,  10  and  24  degrees. 

Prediction  of  the  irpsure  recovery  distributions  using  the  UIM 
method  [2]  will  be  accomp l i si.ed  shortly. 

Measurement  Techniques.  Both  channels  of  pulsed  wire  anemometry  and 
probes  were  delivered,  and  preliminary  experiments  to  check  their  oper- 
ating characteristics  have  commenced.  These  instruments  were  purchased 
at  no  cost  to  the  project. 

A wall-flow  direction  indicator  has  been  developed  in  our  labora- 
tory. The  probe  holds  three  fine  wires,  0.25  inches  long,  that  are 
mounted  parallel  to  each  other  at  a spacing  of  0.1  inch.  The  probe 
plug  fits  flush  into  the  wall  of  a flow  channel  so  that  the  wires  are  held 
perpendicular  to  the  through  flow  direction.  The  two  outer  wires  are 
set  out  from  the  wall  about  1 mm  on  short  prongs  and  the  center  wire 
is  set  in  a shallow  groove.  The  center  wire  is  driven  with  a steady 
current  so  that  it  produces  a hot  wake  which  is  swept  downstream,  or 
upstream,  over  one  of  the  outer  wires.  The  outer  wires,  operated  as 
resistance  thermometers,  detect  the  wake  and,  by  use  of  a simple  com- 
parator circuit,  give  the  instantaneous  sign  of  the  fluid  velocity  in  a 
thin  layer  near  the  wall.  The  system  has  been  successfully  tested  in  a 
reattaching  turbulent  shear  flow,  and  will  be  applied  for  the  detection 
of  the  unsteady  flow  separation  points  in  our  diffuser  experiments. 
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Introduction 


The  overall  goal  of  this  research  program  is  to  provide  a better  under- 
standing of  stall-related  phenomena  in  axial-flow  compressors.  The  aspects 
of  compressor  performance  that  are  of  interest  include  the  onset  of  stall, 
loss  in  performance,  and  the  flow  instabilities  associated  with  stall. 

The  program  involves  both  experimental  and  analytical  efforts.  A 
primary  feature  of  the  experimental  work  is  the  measurement  of  pressures 
directly  on  the  rotor  of  test  compressors,  for  flow  conditions  up  to  and 
including  stall.  For  high-frequency-response  measurements,  special  radio 
telemetry  data  transmission  equipment  has  been  developed  for  use  with  blade- 
mounted  transducers.  Average  and  slowly- varying  pressure  measurements 
are  made  employing  a pressure  scanner  that  has  been  adapted  to  rotate  with 
the  compressor  rotor.  Ports  on  the  compressor  blades  are  connected  by  tubing 
to  the  scanning  valve.  Pressure  measurements  are  made  by  a single  transducer 
located  at  the  shaft  center. 

Experimental  work  is  conducted  employing  a low-speed,  single  stage 
research  compressor  operating  at  approximately  2400  rpm,  and  a recently- 
completed  high-speed  drive  facility.  A three-stage  research  compressor 
designed  for  operation  in  the  13000-17000  rpm  range  is  presently  being 
readied  for  use  in  this  facility. 

During  the  present  reporting  period,  the  use  of  high  response  dynamic 
pressure  probes  to  supplement  on-rotor  pressure  measurements  was  investigated. 
Two  probes  were  installed  downstream  of  the  rotor  of  the  low-speed  test 
compressor,  and  measurements  were  made  with  rotating  stall  present.  Work 
is  proceeding  to  install  similar  probes  in  the  first-stage  stator  row  of 
the  high-speed  test  compressor.  This  test  compressor  is  being  prepared  for 
a series  of  experiments  to  determine  the  stall  point  of  the  rotor  blades, 
and  to  compare  the  observed  behavior  with  predictions  of  a separated-flow 
theory  which  has  been  devised.  Progress  was  made  in  the  development  of  this 
theory,  and  in  several  areas  related  to  the  initial  operation  of  the  high- 
speed research  compressor. 
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Discussion 


Development  of  techniques  for  on-rotor  pressure  measurements  in  low- 
speed  axial-flow  compressors  has  been  previously  reported  in  connection 
with  the  present  program.  In  order  to  study  the  time-response  of  the  rotor 
blades  to  disturbances  which  produce  stall,  it  is  necessary  to  provide  for 
high-frequency-response  pressure  measurements  in  the  stationary  components 
of  the  compressor  as  well. 

Two  dynamic  pressure  probes  (KULITE  XCQH-152-15D)  were  installed  in 
the  single-stage  low-speed  research  compressor  at  120°  circumferential 
separation,  and  on  traverses  which  permitted  radial  movement.  The  probes 
were  located  approximately  15  mm  behind  the  rotor  of  the  compressor.  The 
probes  proved  very  effective  for  indicating  the  presence  of  rotating  stall, 
and  for  measuring  the  radial  and  circumferential  extent  of  the  cells  [1]. 
Preparations  are  underway  for  similar  high-response  flow  measurements  down- 
stream of  the  first-stage  rotor  of  the  high  speed  research  compressor. 

Primary  considerations  involve  the  method  of  mounting  and  radially  traversing 
the  probes,  and  the  flow  blockage  effects  of  the  probes  when  mounted  within 
the  first-stage  stator  now.  The  high-response  probes  will  be  used  to  determine 
the  type  of  unstable  behavior  which  occurs  in  the  rotor,  and  to  measure  the 
rate  of  flow  fluctuations.  A related  experimental  effort  planned  for  the 
high-speed  facility  is  directed  at  determining  the  mean  flow  performance  of 
the  research  compressor  up  to  stall,  for  comparison  with  the  basic  analytical 
model . 

Recent  progress  in  both  efforts  includes  completion  of  a new  test  rotor, 
final  design  of  the  inlet  and  exit  ducting,  design  of  the  first  experimental 
program,  and  selection  of  the  necessary  instrumentation.  The  test  compressor 
has  been  prepared  in  a three-stage  configuration  with  new  bearings,  balanced, 
and  is  now  ready  for  installation  in  the  facility.  The  inlet  and  exit 
ducting,  which  includes  a variable  volume  discharge  plenum  and  quick-release 
valve,  has  been  designed  and  is  presently  under  construction.  The  first 
experiments  will  include  flow  measurements  at  the  compressor  inlet,  after 
the  first-stage  rotor,  after  the  last  (third)  stage,  and  in  the  plenum.  An 
automatic  traversing  mechanism  will-  position  the  measurement  probe  behind 
the  first-stage  rotor. 

The  analytical  effort  is  directed  toward  developing  a basic  model  for 
predicting  stall  in  axial-flow  compressors  which  can  be  compared  with  the 
experiments.  Recent  progress  in  this  aspect  of  the  program  includes  a basic 
procedure  for  calculating  separated  boundary  layers,  a prediction  method  for 
turbulent  boundary  layers  with  separated  flow,  a computer  program  for  radial 
equilibrium,  and  a blade-to-blade  program  for  the  compressible,  inviscid 
flow.  The  procedure  for  calculating  separated  flow  [2]  involves  a simultaneous 
solution  of  the  boundary  layer  and  freestream  at  each  longitudinal  position. 

The  solution  begins  at  the  upstream  boundary  and  iterates  over  the  flow 
field  to  satisfy  the  elliptic  nature  of  the  freestream  equation.  Stability 
and  convergence  of  the  basic  procedure  were  first  demonstrated  with  laminar, 
internal  flow. 

The  turbulent  boundary  layer  prediction  method  is  an  integral  method 
that  can  be  easily  used  with  the  above  separated  flow  procedure.  The  method 
employs  the  momentum  and  kinetic  energy  integral  equations  and,  except  for  the 


22 


extension  to  separated  flow,  is  not  greatly  different  from  a number  of 
previous  methods.  The  velocity  profiles  are  made  up  of  the  usual  logarithmic 
and  wake  functions,  but  in  the  separated  region  the  wake  function  does  not 
extend  to  the  wall.  Thus  in  the  limit,  the  profiles  result  in  a free  shear 
layer  at  a large  distance  from  the  wall,  similar  to  the  Stewartson  profiles 
for  laminar  flow.  The  dissipation  integral  is  being  determined  from  available 
experimental  results,  which  are  quite  limited  for  separated  flow. 

The  radial  equilibrium  program  is  a straight-forward  calculation  that 
can  be  used  with  the  blade- to-blade  program  to  calculate  the  full  blade 
passage.  Streamline  shift  is  determined  from  simple  radial  equilibrium, 
and  the  axial  shape  of  the  stream  sheets,  which  are  assumed  axially 
symmetric,  is  approximated  by  Horlock's  experimental  function. 

The  bl ade-to-blade  program  has  been  completed  for  compressible,  inviscid 
flow  and  is  similar  to  the  Katsanis  program.  However,  the  present  computer 
program,  which  employs  line-overrelaxation,  has  been  developed  specifically 
to  incorporate  a simultaneous  boundary  layer  calculation  that  includes 
separated  flow. 
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Introduction 


The  long-term  objective  of  this  research  program  is  to  ascertain  tlie 
role  which  oncoming  turbulence  plays  in  reducing  the  aerodynamic  losses 
in  flow  through  a blade  cascade  of  an  axial -flow  compressor  at  moderate 
Reynolds  numbers  of  order  of  2x10'’  or  smaller.  At  these  Reynolds  numbers 
[jroliibit ivcly  high  losses  and  even  fully  stalled  blades  arc  induced  by 
laminar  separation  of  the  profile  boundary  layer.  Supjily  of  oncoming  tur- 
bulence of  sufficient  energy  concentrated  at  scales  commensurate  with  the 
thickness  of  the  prevalent  profile  boundary  layer  can  forestall  the  lami- 
nar separation.  Suitable  management  of  the  turbulent  energy  distribution 
possesses  further  the  potential  to  even  generate  atid  sustain  a fully 
attached  turbulent  boundary  layer  on  the  profile  suction  side.  Accumula- 
tion of  turbulent  energy  at  desired  scales  can  be  produced  by  selective 
amplification  of  turbulence.  This  selective  turbulent  energy  intensifi- 
cation is  governed  by  the  vortex-stretching  mechanism  characteristic  to 
forward  stagnation  flow. 

The  research  program  is  divided  into  three  phases  of  increasing  com- 
plexity for  the  sake  of  ensuring  its  methodical  and  successful  completion. 
In  all  these  three  phases  the  evolution  of  the  oncoming  turbulent  energy, 
its  selective  amplification  and  its  effects  on  the  body  boundary  layer 
are  to  be  investigated.  The  bodies  to  be  utilized  are:  (1)  a circular 
cylinder  in  the  first  phase;  (2)  an  isolated  airfoil  in  the  second  phase; 
and,  (3)  a stationary  blade  cascade  in  the  last  phase.  Essentially,  the 
first  phase  represents  a diagnostic  study  regarding  the  flow  features  of 
interest  and  the  foundations  for  putting  forward  an  adequate  theoretical 
model . 


Discussion 


The  current  research  efforts  concentrate  on  investigating  the  oncom- 
ing turbulence  evolution  and  the  effects  of  the  amplified  turbulence  upon 
the  flow  about  a circular  cylinder,  i.e.,  the  first  phase  of  the  research 
program.  An  extensive  visualization  study  of  the  flow  near  the  stagnation 
zone  of  a circular  cylinder  was  conducted  and  completed.  The  results  of 
this  visualization  investigation  are  reported  in  a Project  SQUID  technical 
report  [l]  which  was  published  recently.  A summary  movie  titled  "Vorti- 
city  Amplification  in  Stagnation  Flow"  which  incorporates  the  most  instruc- 
tive views  of  the  flow  was  further  produced. 

The  visualization  supplied  an  in-depth  qualitative  apperception  of 
the  flow.  A frame  by  frame  scrutiny  of  selected  movie  strips  led  further 
to  the  acquisition  of  a reasonable  quantitative  interpretation  of  the 
gross  flow  structure.  Kssentially , the  visualization  study  provided  sig- 
nificant and  explicit  evidence  concerning:  (1)  the  selective  stretching 
of  cross-vortex  tubes;  (2)  the  streamwise  biased  tilting  of  stretched 
cross-vortex  tubes;  (3)  the  existence  of  a coherent  vortex  flow  structure 
near  the  stagnation  zone;  (4)  the  interaction  of  the  amplified  turbulence 
with  the  cylinder  laminar  boundary  layer;  and,  (5)  the  fostering  and 
growth  of  a turbulent  boundary  layer.  The  visualization  revealed  clearly 
that  the  cross-vortex  tubes  conveyed  by  the  diverging  stagnation  flow  con- 
stitute an  organized  substructure  within  the  overall  turbulent  field 
which  is  triggered  to  its  fullest  manifestation  by  the  stretching  mecha- 
nism. 


The  effect  of  the  amplified  turbulence  upon  the  position  of  the  sepa- 
ration line  was  further  visualized.  By  and  large,  it  was  found  that  the 
separation  line  can  be  adjusted  within  a range  of  7 to  beyond  10  cm  (3 
to  4 in)  along  tlic  cylinder  circumference  depending  upon  the  character- 
istics of  the  oncoming  turbulence.  A short  report  describing  the  results 
of  this  endeavor  is  currently  l)eing  prepared. 

A first  survey  of  the  turbulence  evolution  along  the  stagnation 
streamline,  i.e.,  along  the  X2-axis,  was  carried  out  by  means  of  a 
single  hot-wire  anemometer.  Measurements  were  performed  at  six  Reynolds 
numbers,  viz.,  at  Re^  = 5x10**,  7x10**,  9x10**,  1.2x10^,  1.6x10^  and  2x10^, 
whore  Rep  is  the  cylinder  Reynolds  number  based  on  the  freostream  velo- 
city and  the  cylinder  diameter.  At  each  Reynolds  number  the  measurements 
were  conducted  at  21  stations  ranging  From  14  radii  upstream  of  the  cylin- 
der up  to  0.01  radii  close  to  the  cylinder.  The  surveys  were  carried  out 
with  a turbulence-generating  grid  made  up  of  vertical  rods  [l]  installed 
14.64  radii  upwind  of  the  cylinder.  During  this  first  survey  the  stream- 
wise  turbulent  velocity  (rms) , i.e.,  the  U2-component , and  its  power  spec- 
tral density  function  were  measured  on  line.  In  addition,  the  data  was 
recorded  on  I'M  magnetic  tape.  The  reduction  and  analysis  of  the  amassed 
data  is  currently  being  conducted.  A preliminary  analysis  of  the  data 
clearly  revealed  that  the  turbulent  energy  undergoes  manifold  amplifica- 
tion close  to  the  l)ody  at  scales  larger  than  the  neutral  one.  Amplifica- 
tion ratios  greater  than  even  100  times  were  obtained  near  the  outer  edge 
of  the  tl)eoretical  laminar  boundary  layer. 
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The  construction  of  a single  NACA  65-010  airfoil  was  completed.  This 
airfoil  has  a chord  of  122  cm  (48  in).  Along  its  circumference  120  pres- 
sure taps  were  drilled.  This  airfoil  will  be  utilized  during  the  second 
phase  of  this  research  program. 

I’urther  efforts  regarding  tlic  matching  of  the  inner  and  outer  solu- 
tions of  the  vorticity  amplification  theory  have  been  pursued.  Numerical 
schemes  for  carrying  out  the  matching  by  means  of  a composite  asymptotic 
expansion  are  being  evaluated.  An  analysis  of  the  differences  and  the 
similarities  between  the  rapid  distortion  theory  and  the  vorticity  ampli- 
fication theory  has  been  initiated.  The  main  objective  of  this  analysis 
is  the  ascertainment  of  the  effect  of  the  turbulent  velocity  gradient  upon 
the  stretching  of  cross-vortex  tubes.  An  additional  aim  is  the  examina- 
tion of  the  role  of  viscous  dissipation  with  regard  to  the  stretching 
mechanism. 
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Introduction 


The  problem  of  turbulent  boundary  layer  separation  due  to  an 
adverse  pressure  gradient  is  an  important  factor  in  the  design  of 
many  devices  such  as  jet  engines,  rocket  nozzles,  airfoils  and 
helicopter  blades,  and  the  design  of  fluidic  logic  systems.  Until 
the  last  three  years  little  quantitative  experimental  information 
was  available  on  the  flow  structure  downstream  o^"  separation  be- 
cause of  the  lack  of  proper  instrumentation. 

In  1974  after  several  years  of  development,  a one  velocity 
component  directionally-sensitive  laser  anemometer  system  was 
used  to  reveal  some  new  features  of  a separating  turbulent  bound- 
ary layer  [1].  The  directional  sensitivity  of  the  laser  anemo- 
meter system  was  necessary  since  the  magnitude  and  direction  of 
the  flow  must  be  known  when  the  flow  moves  in  different  directions 
at  different  instants  in  time  [2].  In  addition  to  much  turbulence 
structure  information,  it  was  determined  (1)  that  the  law-of-the- 
wall  velocity  profile  is  apparently  valid  up  to  the  beginning  of 
intermittent  separation;  (2)  that  the  location  of  the  beginning  of 
intermittent  separation  or  the  upstreammost  location  where  separa- 
tion occurs  intermittently  is  located  close  to  where  the  free- 
stream  pressure  gradient  begins  to  rapidly  decrease;  (3)  that  the 
normal  stress  terms  of  the  momentum  and  turbulence  kinetic  energy 
equations  are  important  near  separation;  and  (4)  that  the  separated 
flowfield  shows  some  similarity  of  the  streamwise  velocity  U,  of 
the  velocity  fluctuation  u',  and  of  the  fraction  of  time  that  the 
flow  moves  downstream  [3]. 

Based  upon  these  results,  modifications  [4,5]  to  the  Bradshaw, 
et  al . [6]  boundary  layer  prediction  method  were  made  with  signi- 
ficant improvements.  However,  this  prediction  effort  pointed  to 
the  need  to  understand  the  relationship  between  the  pressure  gra- 
dient relaxation  and  the  intermittent  separation  region  structure. 
Another  limiting  factor  for  further  refinement  of  the  prediction 
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of  separated  flows  is  the  lack  of  fundamental  velocity  and  turbu- 
lence structure  information,  especially  in  the  backflow  region. 

Thus,  the  objective  of  the  current  research  program  is  to  provide 
this  information  by  using  a directionally-sensitive  laser  anemo- 
meter system  to  determine  quantitatively  the  turbulence  structure 
of  a separating,  separated,  and  reattached  turbulent  boundary  layer. 


Discussion 


This  current  research  program  was  begun  October  1,  1976,  to 
obtain  laser  anemometer  measurements  of  the  separating  flow  of  another 
adverse  pressure  gradient  turbulent  boundary  layer  for  an  airfoil 
or  cascade  blade  type  pressure  distribution.  As  stated  in  the  last 
report,  considerable  effort  has  been  made  to  avoid  mean  flow  three- 
dimensionality.  Specially  designed  wall  suction  and  tangential 
wall  jet  boundary  layer  controls  and  peripheral  equipment  were  in- 
stalled into  the  wind  tunnel  test  section.  The  flow  produced  by 
these  controls  was  determined  to  be  two-dimensional  within  1%.  Recent 
measurements  indicate  that  the  momentum  from  the  specially-designed 
wall  jets  persists  farther  dc'^nstream  than  other  more  conventional 
designs.  Direct  measurements  of  U and  V in  the  backflow  region  with 
the  directly-sensitive  laser  anemometer  system  have  proved  that  the 
time-averaged  two-dimensional  continuity  equation  is  satisfied. 

A number  of  velocity  profiles  have  been  measured  upstream  of 
separation  with  a hot-wire  anemometer.  The  streamwise  skin-friction 
behavior  is  well  described  by  the  Ludwieg-Tillmann  relation.  The 
mean  and  streamwise  fluctuation  velocity  profiles  agree  with  pre- 
viously published  results. 

Cin^  films  of  laser  illuminated  smoke  have  clearly  revealed 
that  the  large  eddy  structure  supplies  most  of  the  near  wall  back- 

flow.  Many  measurements  of  U,  V,  u^,  v^,  (U  - V )//?”,  -uv  the 
flow  reversal  intermittency,  the  skewness  and  flatness  of  the  velocity 
probability  distributions,  and  velocity  spectra  are  being  obtained 
in  this  region  with  the  laser  anemometer  system.  The  signal  quali- 
ty of  this  backscattering  system  has  been  considerably  upgraded 
in  the  last  several  months.  Up  to  10^  data  samples  per  second  are 
being  obtained. 

Three  distinct  layers  in  the  backflow  have  been  observed  so 
far.  Nearest  the  wall  is  a viscous  layer  that  is  driven  by  the 

large-scaled  low  frequency  unsteadiness  of  u^  and  v^  of  the  outer 

flow.  In  this  layer,  is  about  twice  |U|.  The  Reynolds  shearing 

stress  -uV  is  immeasurably  small,  indicating  that  the  mean  velocity  of 
this  viscous  layer  is  not  governed  by  turbulent  momentum  transport. 

It  appears  that  this  layer  behaves  as  an  unsteady  viscous  boundary 
layer  that  is  driven  by  u and  v fluctuations  in  the  second  layer. 
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The  normal  stresses  term  a (u^  - v^)^x  in  the  momentum  equation 

and  (u^  - v^)  ^L)^x  in  the  energy  equation  [4]  are  not  insignificant 
in  the  both  the  viscous  layer  nearest  the  wall  and  in  the  second  layer. 
The  mean  velocity  profile  for  U in  the  second  layer  is  amazingly  flat, 

i.e.,  aU/Sy  = 0,  which  is  in  substantial  agreement  with  the  near  wall 

measurements  of  Simpson  et  al . [3].  -uv  is  also  immeasurably  small 
in  this  layer.  It  is  not  surprising  then  that  the  mean  velocity  profile 
of  the  viscous  layer  nearest  the  wall  looks  much  like  a Blasius  flat 
plate  laminar  layer. 

The  third  and  outermost  layer  of  the  mean  backflow  is  really 
part  of  the  forward  flow  of  the  outer  region.  The  Reynolds  stresses 
are  growing  and  the  mean  velocity  profile  is  taking  on  a more  mixing 
layer  like  behavior. 

Further  measurements  are  being  made  and  analytical  explanation 
of  this  behavior  is  being  pursued. 
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Introduction 


The  study  of  oxidation  reactions  in  shock  tubes  has  been  stimulated 
by  the  use  of  fast,  accurate  numerical  integration  routines.  Now  it  is 
possible  for  kineticists  to  more  definitively  test  various  oxidation  mechan- 
isms by  a detailed  comparison  of  calculated  and  observed  concentration-time 
profiles.  Although  the  H2/02/Ar  and  H2/02/C0/Ar  systems  have  been  success- 
fully studied  by  this  approach,  extension  to  even  simple  hydrocarbon  systems 
like  CH4/O2  has  been  limited  by  lack  of  reliable  high  temperature  rate  con- 
stants. A common  practice  has  been  to  extrapolate  low  temperature  flow  sys- 
tem data  to  the  temperature  range  of  interest.  Unfortunately  this  approach 
can  lead  to  serious  errors;  recent  studies  have  convincingly  demonstrated 
that  many  reactions  of  importance  in  combustion  mechanisms  exhibit  markedly 
"non-Arrhenius"  rate  constants.  In  this  light  it  appears  to  be  most  de- 
sirable to  measure  rate  constants  in  the  same  high  temperature  regime  where 
they  will  be  used  to  test  the  combustion  mechanisms.  However,  it  is  equally 
important  that  these  data  be  obtained  from  relatively  simple  systems  where 
assignment  of  the  desired  rate  constant  is  not  contingent  upon  proper  assign- 
ment of  a complex  mechanism  and  the  associated  rate  constants. 

One  such  system  results  from  the  substitution  of  N2O  for  O2  in  combustion 
studies.  Recent  work  in  this  laboratory  [1]  showed  that  N2O  is  a particularly 
useful  source  of  oxygen  atoms  between  2000-3000  K.  Thus,  a study  of  combustion 
systems  where  N2O  replaced  O2  should  provide  useful  information  about  rates  of 
oxygen  atom  reactions  at  high  temperatures.  The  primary  advantage  of  N2O  as 
an  oxidant  is  that  oxygen  atom  reactions  will  occur  in  an  environment  where 
the  concentration  of  molecular  oxygen  is  much  less  than  a normal  combustion 
system;  this  considerably  simplifies  the  kinetic  analysis.  Prudence  dictates 
that  such  a substitution  first  be  tested  on  a known  system.  For  this  reason 
our  first  efforts  in  the  SQUID  program  utilized  hydrogen  as  the  fuel  molecule. 
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The  hydrogen  work  has  been  completed  [2]  and  suggests  the  substitution  of 
N2O  for  O2  is  a viable  method  to  measure  rates  of  oxygen  atom  reactions  at 
high  temperatures.  Furthermore,  it  appeared  that  the  postulated  mechanism 
for  the  N2O  system  was  reasonable.  Data  were  then  collected  on  a variety  of 
N2O/CH4  mixtures.  Analysis  of  these  data  at  high  temperatures  (T  > 2400  K) 
yielded  values  of  the  rate  constant  for  0 + CH4  ->■  CH3  + OH  in  good  agreement 
with  recently  measured  high  temperature  values.  It  is  significant  to  note 
that  these  values  are  over  five  times  larger  than  those  obtained  by  extrapo- 
lation of  low  temperature  data;  this  observation  reinforces  the  thesis  that 
reaction  rates  used  in  combustion  mechanisms  should  come  from  high  tempera- 
ture measurements.  Detailed  comparisons  of  calculated  and  observed  profiles 
at  temperatures  less  than  2200  K indicated  unsuspected  complexities  in  the 
N2O/CH4  mechanism.  Subsequent  experiments  suggested  that  reactions  of  formal- 
dehyde (CH2O)  were  responsible.  Thus  our  efforts  have  shifted  to  a study  of 
these  reactions. 


Discussion 


During  the  last  six  months,  a variety  of  CH20/Ar  and  CH2O/N2O  mixtures 
have  been  investigated.  Particular  emphasis  has  been  placed  upon  the  CH20/Ar 
system  [3].  Here  CH2O  decay  was  monitored  in  a series  of  mixtures  where  the 
formaldehyde  concentration  varied  from  0.1%  to  1.0%.  These  data  suggest  that 
the  reaction  is  complex,  with  secondary  reactions  playing  a significant  role 
in  the  overall  disappearanceof  CH2O.  A possible  mechanism  is  the  following: 


CH2O  + M -V  HCO  + H + M (1) 
CH2O  + H HCO  + H2  (2) 
HCO  + M H + CO  + M (3) 
HCO  + H H2  + CO  (4) 


At  present,  we  are  trying  to  quantitatively  model  the  decay  behavior,  as 
a function  of  both  temperature  and  concentration,  in  terms  of  this  mechanism. 

Experiments  have  also  been  performed  in  which  the  formaldehyde  decay  is 
measured  in  CH2O/N2O  mixtures.  As  expected,  the  observed  decay  rate  is  en- 
hanced when  N2O  is  added.  The  increase  is  probably  due  to  the  fact  that  the 
additional  reaction  channels 


0 + CH2O  ->  HCO  + OH 
OH  + CH2O  > HCO  + H2O 

are  now  available.  Successful  analysis  of  this  system  (coupled  with  earlier 
work  on  0-atom  production  and  CO2  production  in  CH20/N20/C0/Ar  mixtures) 
should  yield  rate  constants  for  these  reactions. 
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When  these  formaldehyde  reactions  are  characterized,  we  should  then  be  able 
to  proceed  rather  quickly  in  the  CH4/N2O  analysis.  Of  particular  interest  here 
is  that  the  formaldehyde  work  should  be  applicable  to  both  the  CH4/N2O  and 
CH4/O2  systems.  Thus,  the  work  on  CH4/N2O  should  yield  detailed  mechanistic 
information  directly  applicable  to  practical  combustor  systems  as  well  as  the 
desired  0 + CH4  ->■  CH3  + OH  rate  constant. 


Notes  and  References 


1.  A.M.  Dean  and  D.C.  Steiner,  J.  Chem.  Phys.,  76^,  598  (1977).  (Project 
SQUID  Technical  Report  UMO-l-PU.  Earlier  N2O  work  is  referenced  here.) 

2.  A.M.  Dean,  D.C.  Steiner,  and  E.E.  Wang,  Project  SQUID  Technical  Report 
UM0-2-PU.  (To  be  published  in  Combustion  & Flame.) 

3.  A.M.  Dean  and  E.E.  Wang,  "Shock  Tube  Studies  of  Formaldehyde  Pyrolysis" 
(Extended  abstract  accepted  for  presentation  at  the  Seventeenth  Combustion 
Symposium) . 
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COMBUSTION  KINETICS  AND  REACTIVE  SCATTERING  EXPERIMENTS 


Yale  University,  New  Haven,  Connecticut 
Subcontract  No.  4965-16 


J.  B.  Fenn,  Principal  Investigator 
B.  Halpern  and  M.  Labowsky 

Introduction 

The  combustion  of  hydrocarbon  fuels  has  been  man's  most  used  source 
of  useful  energy  for  much  of  this  century.  The  chemical  reactions  which 
it  involves  have  been  among  the  most  studied.  And  yet,  there  remains  un- 
certainty as  to  the  nature  of  the  first  reactive  step  in  the  complex 
sequence  of  reactions  by  which  oxygen  and  hydrocarbon  molecules  become  hot 
combustion  products.  This  investigation  comprises  an  attempt  to  identify 
that  first  reactive  event  and  to  determine  its  cross  section  by  means  of 
molecular  beam  scattering  methods.  The  prospective  advantage  of  such 
methods  is  that  they  can  examine  the  consequences  of  a single  collision 
between  individual  molecules.  By  the  same  token  they  are  substantially 
limited  in  their  ability  to  probe  intermediate  reaction  steps  which  involve 
species  of  transient  existence  such  as  free  radicals  not  readily  obtainable 
as  beams.  In  addition  to  this  new  venture  in  combustion  kinetics  we  have 
been  continuing  a study  of  the  evaporation  and  combuocion  of  arrays  of 
droplets.  This  study  is  based  on  an  adaptation  of  the  method  of  images 
which  has  been  successful  in  solving  Laplace's  equation  as  it  applies  to 
electrostatic  problems  involving  arrays  of  charged  particles. 


A.  Reactive  Scattering.  The  objective  is  to  identify  and  measure 
cross  sections  for  the  first  reactive  steps  in  the  combustion  oxidation 


of  paraffin  hydrocarbons.  In  order  to  avoid  ambiguity  we  are  trying  to 
carry  out  molecular-beam-like  experiments  so  that  we  can  be  sure  that  any 
product  we  see  is  the  result  of  a single  collision.  Our  approach  is  to 
employ  uncollimated  opposing  beams  of  oxygen  and  hydrocarbon  molecules 
comprising  free  jets  from  small  sonic  nozzles  exhausting  into  an  evacuated 
region.  After  a collision  both  reactant  and  product  molecules  are  collected 
on  a cryogenic  trap.  Collection  can  continue  for  suitably  long  periods  of 
time  so  that  detectable  amounts  of  product  will  be  obtained  even  if  the 
reaction  cross  section  is  very  small.  The  reaction  chamber  is  then  isolated 
and  heated  so  that  the  trapped  species  return  to  the  gas  phase  and  can  be 
swept  out  by  a stream  of  helium  for  analysis  by  gas  chromatography. 

We  have  built  the  reaction  system  and  have  tried  a number  of  runs 
without  yet  obtaining  evidence  of  reaction.  The  absence  of  product  may  mean 
that  the  cross  section  is  too  small  for  us  to  detect.  We  think  that  the 
likely  product  may  be  the  olefin  corresponding  to  the  paraffin  feed,  i.e., 
butene  from  butane.  Because  the  separation  of  these  two  peaks  is  not  all 
that  great,  the  apparent  absence  of  any  product  may  be  due  to  the  "loss" 
of  a small  product  butene  peak  in  the  large  reactant  butane  peak.  Con- 
sequently we  have  been  trying  to  carry  out  partial  fractionation  of  the 
reactant-product  mixture  before  analysis.  So  far  we  have  been  unsuccessful. 
Therefore,  we  are  now  preparing  to  run  a reaction  which  should  have  a larger 
cross  section  and  whose  products  should  be  more  easily  resolved  on  the 
chromatograph.  Our  candidate  is  the  "oxidation"  of  hydrocarbon  by  chlorine 
to  form  HCl  and  an  alkyl  chloride. 

B.  Evaporation  and  Combustion  of  Droplet  Arrays.  In  previous  work  we 
have  used  the  method-of-images  to  solve  Laplace's  equation  for  the  combustion 
of  arrays  of  interacting  but  slowly  evaporating  fuel  particles.  The  results 
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indicated  that  under  typical  conditions  for  spray  combustion  the  droplets 
burn  in  a group  mode  with  a common  flame  rather  than  as  aggregations  of 
single  particles  each  with  its  own  flamelet.^^^  It  became  of  interest, 
therefore,  to  extend  the  method  to  the  case  of  rapidly  evaporating  droplets. 
In  this  case,  the  convective  terms  associated  with  vapor  transfer  (Stefan 
flow)  cannot  be  neglected.  The  presence  of  these  terms  makes  the  rapid 
evaporation  transport  equations  non-linear.  Hence,  they  cannot  be  solved 
directly  by  superposition. 

By  means  of  suitable  transformations  we  have  reduced  the  equations  to 

a form  which  can  be  solved  by  the  Images  method  as  previously  described. 

It  turns  out  that  the  particle-particle  interactions  can  substantially 

reduce  the  burning  rates  of  fuel  droplets  even  for  fairly  large  particle 

spacings,  i.e.,  of  the  order  of  30  particle  radii.  Moreover,  the  so-called 
2 

"D  -Law",  which  asserts  that  the  surface  area  of  an  isolated  burning  droplet 
decreases  linearly  with  time,  does  not  hold  in  tlie  case  of  Interacting 
droplets. 

References 

1.  M.  Labowsky  and  I).  E.  Rosner,  Proceedings  of  the  Symposium  on  Evaporation 
and  Combustion  of  Fuel  Droplets,  American  Chemical  Society  Centennial 
Meeting,  San  Francisco  (1976),  in  ACS  Advances  in  Chemistry  (in  press). 
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No  work  has  been  undertaken  in  this  project  since  October 
in  view  of  funds  for  extension  of  contract  not  being  in  hand. 
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III.  MEASUREMENTS 


TURBULENT  STRUCTURE  DETERMINATION  BY  RAMANOGRAPHY 


Yale  University,  New  Haven,  Connecticut 
Subcontract  No.  8960-29 


Professor  R.  K.  Chang,  Principal  Investigator 
Mr.  M.  Long,  Assistant  in  Research 
Mr.  B.  Webber,  Assistant  in  Research 


Introduction 


The  purpose  of  our  research  is  to  develop  new  spectroscopic 
techniques  that  are  capable  of  providing  fundamental  information 
in  the  areas  of  turbulence  and  combustion.  Presently,  we  are  inves- 
tigating approaches  to  measure  the  instantaneous  concentration  pro- 
file of  specific  gases  in  a jet,  the  instantaneous  temperature  within 
a flame,  and  the  velocity  flow  field  in  a jet.  All  three  of  these 
optical  techniques  are  being  pursued  on  a well  calibrated  nozzle  and 
a steady  non-luminous  flame. 


Discussion 


Concentration  Profile.  The  concentration  profile  of  a specific  gas 
is  determined  by  measuring  the  Raman  scattering  intensity  of  that 
gas  within  a one-dimensional  track  defined  by  the  laser  radiation. 
Instantaneous  mapping  of  the  Raman  intensity  along  this  track  is 
accomplished  by  using  a pulsed  laser  to  provide  the  high  energy  in- 
cident radiation  and  a low-light  level  TV  camera  to  detect  the  weak 
Raman  scattered  radiation  which  has  been  isolated  from  the  elasti- 
cally scattered  radiation  by  a double  monochromator  or  an  inter- 


ference filter  centered  at  the  Raman  wavelength.  In  order  to  freeze 
any  small  scale  turbulent  structure  (10~^  cm)  with  a flow  field 
(velocity  between  10^  - IC^  cm/sec) , the  laser  pulse  duration  needs 


to  be  as  short  as  10  sec.  However,  the  requirement  for  the  mini- 
mum number  of  Raman  photons  for  a reasonable  signal-to-nolse  ratio 
demands  that  the  ruby  laser  has  about  10  Joules  per  pulse.  We  are 
developing  a nonlinear-optical  method  to  stretch  a conventional  Q- 
switched  ruby  pulse  from  50  x 19“^  sec  to  10“”  sec.  Once  the  Raman 
intensity  profile  is  in  the  computer,  we  can  calculate  its  auto- 
correlation, which  should  be  Identical  to  the  density  autocorrela- 
tion function  within  a time  Interval  defined  by  the  pulsed  laser 
duration . 

Temperature  Measurement.  The  temperature  averaged  along  a laser 
track  passing  through  a flame  is  determined  by  detecting  the  spec- 
tral llneshape  of  the  Stokes  Raman  or  anti-Stokes  Raman  radiation 
emanating  along  this  track.  The  Stokes  Raman  lineshapes  consist 
of  a series  of  peaks  arising  from  vibration  transitions  V=0’-1, 

V=l'-2,  and  V=2'-3,  while  the  anti-Stokes  Raman  lineshapes  corres- 
pondingly consist  of  V=1>0,  V=2'-l,  and  V=3>2  peaks.  Kurthermore, 
the.  asymmetrical  1 ineshape  of  each  peak  contains  information  on  the 
rotational  temperature.  Instantaneous  determination  of  the  tempera- 
ture within  an  Il2-air  porous  plug  liurner  is  achieved  by  using  a 
pulsed  laser  as  the  incident  radiation  source,  the  spectrograph  to 
disperse  tlie  spectral  content  of  the  Raman  radiation  (Stokes  or 
anti-Stokes),  and  a low-light  level  TV  camera  to  record  the  entire 
Raman  spectral  1 ineshape  from  the  laser  track  which  passes  through 
the  flame.  The  viliratlonal  and  rotational  temperature  profiles 
along  this  track  can  he  deduced  by  curve  fitting  the  observed  line- 
shapes  with  the  aid  of  the  computer. 

Velocity  Flow  Field.  The  velocity  flow  field  within  a jet  is  accom- 
plished bv  seeding  the  jet  with  0.3p  particles  and  observing  their 
bright  trajectories  upon  illumination  of  a properly  shaped  light 
pulse  or  a series  of  short  light  pulses.  This  technique  is  a Lag- 
rangian  approach  since  information  is  obtained  by  following  the  tra- 
jectory of  a particle.  The  hot-wire  anemometer  or  the  LDV  is  a 
Eulerian  approacti  where  the  particle  enters  a volume  element  fixed 
in  the  laboratory  frame. 

Immt'diate  Plans.  The  one-d irnensional  concentration  profile  measure- 
ments have  lieen  Initiated  using  a ew  argon  laser  and  a round  nozzle 
emitting  I’reon-12  gas  in  air.  Velocity  flow  field  measurements 
using  tlie  two-dimensional  capaliility  of  our  TV  camera  liavc  also 
lieen  started.  We  slial  1 continue  tlie  concentration  and  velocity 
flow  field  measurements  and  then  combine  this  information  witli  the 
time-resolved  velocity  data  deduced  from  a hot-wire  anemometer  in 
order  to  fully  characterize  the  nozzle  under  various  initial  condi- 
tions . 
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Electronic  conversion  of  our  detection  system  to  record  two- 
dimensional  information  has  been  completed.  This  will  enable  us  to 
map  out  the  two-dimensional  concentration  profile  witliln  a plane 
defined  by  the  laser  beam,  as  well  as  the  one-dimensional  tempera- 
ture profile  along  a laser  track.  Computer  software  needs  to  he 
developed  before  we  can  test  the  sensitivity  of  our  low-light 
level  TV  camera  system,  as  we  know  that  the  Raman  signal  will  be 
reduced  by  the  number  of  resolvable  elements  in  the  second  dimen- 
sion. 
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CARS  INVESTIGATIONS  IN  SOOI’ING  FLAMES 


United  Technologies  Research  Center 
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Alan  C.  Eckbreth,  l^incipal  Investigatfjr 


Introduction 


Jjaser  Raman  spectroscopic  techniques  for  c<jimbustion  diagnostics  have  utide/’gone 
considerable  develojiment  in  the  past  several  years  and  are  now  being  employed  in  a 
variety  of  fundamental  flame  investigations.  Instrumental ly  however,  practical 
combustion  devices  possess  flame  environments  which  differ  markedly  from  those 
tyjiically  (laminar  premixed,  hydrogen  diffusion)  employed  in  fundamental  studies. 
I'ractical  devices  contain  flames  which  can  be  highly  [>articulate  laden  and  hence, 
limiinous,  if  hot,  and  turbulent.  These  conditions  lead  to  a variety  of  severe, 
naturally  occurring  or  laser  induced  interferences  which  must  be  overcome.  Of 
these,  laser  modulated  particulate  incandescence  appears  to  be  the  most  severe. 

When  the  soot  particulate  loadings  become  moderate,  on  the  order  of  10“®  gm/cra^  or 
larger,  the  spontaneous  Raman  signal  to  laser  modulated  soot  incandescence  inter- 
ference ratio  can  become  unacceptably  low  for  measurement  purposes. 

Coherent  anti-Stokes  Raman  spectroscopy  (CARS)  appears  as  an  attractive 
alternative  to  spontaneous  Raman  scattering  for  practical  combustor  diagnosis. 
First,  the  CARS  process  generates  signals  generally  several  orders  of  magnitude 
greater  than  those  possible  with  spontaneous  Raman  scattering.  Second,  the  CARS 
radiation  emerges  as  a coherent  beam  which  can  be  completely  collected  using  high 
f number  collecting  optics  and  spatial  filtering.  This  not  only  leads  to  high 
signal  collection  efficiency,  but  low  interference  collection  as  well  due  to  the 
greatly  reduced  solid  angles  employed.  CARS  signal  strengths  appear  adequately 
large  to  overcome  interferences  from  both  natural  background  luminosity  and  laser 
modulated  particulate  incandescence.  However,  CARS  generation  in  sooting  flames 
had  not  been  demonstrated.  In  sooting  flames,  there  is  the  potential  for  the 
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generation  of  nonlinear  interferences  from  soot  particulates  and  soot  vaporization 
products.  The  objective  of  the  present  research  is  to  address  CARS  generation  in 
sooting  flames,  to  examine  and  develop  techniques  to  suppress  nonlinear  inter- 
ferences, and  to  perform  CARS  temperature  measurements  in  such  flames.  Such  a 
program  logically  precedes  measurement  attempts  in  actual  research  scale  combustors. 


Discussion 


CARS  generation  has  been  studied  in  a laminar,  sooting,  propane  diffusion 
flame.  Such  a flame  was  studied  extensively  earlier  in  conjunction  with  the 
investigations  of  laser  modulated  soot  incandescence  interferences  in  spontaneous 
Raman  diagnostics  (Refs.  1,  2).  The  flame  is  sustained  on  a 0.6  cm  i.d.  stainless 
steel  tube.  The  flame  in  appearance  resembles  that  of  a candle  and  is  highly 
sooting.  iYevious  Mie  measurements  (Ref.  l)  indicated  a soot  number  density  on 
the  order  of  10^*^  cm“3  with  an  average  particle  diameter  of  UooS.  This  is  a very 
high  soot  level  producing  an  attenuate  on  in  transmitted  light  of  a few  tenths  of  a 
percent  f>er  mm  pathlength,  thus  serving  as  a rigorous  test  of  the  capabilities  of 
CARS.  To  obtain  the  requisite  spatial  resolution,  a crossed-beam  phase-matching 
a[)proach  (Ref.  3)  develo])ed  at  UT'RC  is  employed.  The  technique  is  termed  BOXCARS 
based  u[)on  the  sha()e  of  the  phase-matching  diagram,  and  crosses  three  beams,  two 
"pump,"  one  "Stokes"  at  ajjpropriate  angles  to  generate  phase-matched  CARS.  The 
technique  circumvents  the  poor  and  orten  ambiguous  axial  spatial  resolution 
encountered  with  conventional,  collinear  phase -matching.  In  the  initial  experiments, 
the  CARS  was  generated  in  a sample  volume  approximately  0.4  mm  in  dia  by  1 mm  long. 

The  UTRC  CARS  apparatus  employs  a frequency-doubled  neodymium  laser  to  provide  a 
pump  component  at  5320^,  part  of  which  is  split  off  to  generate  a tunable  broadband 
"Stokes"  component  from  a laser  pumped  dye  oscillator/amplifier  combination. 

In  the  first  experiments  performed,  strong  incoherent  and  coherent  inter- 
ferences, roughly  comjiarable  to  the  signal  from  the  ground  vibrational  state  of  N2 
were  encountered.  Incoherent  features  are  generally  suppressible  in  coherent  spec- 
troscopy. Coherent  features  are  obviously  more  problematical.  In  this  instance, 
the  coherent  feature  arose  from  three  wave  mixing  in  laser  produced  C2.  At  the  high 
laser  intensities  used  to  produce  CARS,  significant  soot  vaporization  will  occur. 

This  occurs  on  even  a nanosecond  time  scale  as  a study  of  the  incoherent  emissions 
iiere  demonstrated.  This  is  due  to  an  absence  of  adequate  heat  transfer  processes 
leading  to  a flashed  vajiorj zation  of  the  soot  surface.  A major  pro<luct  of  soot 
vaporization  is  Cm.  CARS  can  l)e  elcctroniealJy  resonantly  enhanced  when  either 
tfie  CARS  ratiiation  itself  or  the  pum(i  component  resides  near  an  electronie  transition. 
Tile  CARD  general, ion  in  the  Laser-  firoduced  Cp  vapor  is  resonantly  enhanced  because 
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of  the  coincidence  of  the  N2  CARS  and  C2  electronic  transition  frequencies.  C2 
has  a major  Swan  band  transition,  aSr^  (v-O)  -»  (v  i),  at  4737*1^  midway 

between  the  No  anti-Stokes  bands  at  4733^  (v  O)  and  4740X  (v  1). 

An  u[)tical  collection  system  was  set  u]!  at  ri^ht  an(/J.cs  to  the  central  axis 
and  the  laser  induced  emissions  examined  thriiUK.h  blockliif'  and  narrowband  inter- 
ference filters.  Both  laser  modulated  incandescences  and  laser  Induced  Swan 
emissions  were  examined.  The  Swan  emissions  were  about  an  order  of  magnitude  more 
intense  than  the  incandescences  and  found  to  be  more  pronounced  with  tlie  broadband 
emission  centered  at  6073^  than  with  532o8  excitation,  probably  due  to  the  presence 
of  three  Swan  absorptions  within  the  broadband  (27O  cm“^  iWIUl)  6073^  Stokes  laser 
band.  This  result  suggested  decreasing  the  dye  laser  spectral  bandwidth  to  decrease 
the  Cp  Swan  absorptions  and  to  increase  the  CARS  signal  strength.  Also  insertion 
of  a polarization  filter  in  the  polarized  CARS  beam  reduces  any  remaining  incoherent 
contributions  in  half.  In  so  doing,  fairly  clean  N2  spectra  were  obtained  which 
exhibited  minimal  Cp  interferences  and  which  permitted  a temperature  measurement 
in  the  flame  of  approximately  2300°K.  The  adiabatic  flame  temperature  for  propane/ 
air  is  2250°K.  The  foregoing  result  is  very  encouraging  since  the  soot  densities 
in  the  propane  flame  are  quite  high  and  is  a favorable  indication  of  the  practical 
viability  of  CARS  diagnostics. 

Kxjieriments  are  continuing  in  this  vein  to  optimize  the  Stokes  laser  bandwidtli 
and  to  mii|)  out  the  temfierature  field  in  a highly  sooting  flame.  The  CARS  teiriperatufe 
measurements  will  most  likely  be  com|)ared  with  those  determined  by  soditun  i.ine  or 
soot  reversal  techniques. 
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Introduction 

Fluctuation  values  of  temperature  have  been  obtained  from  vibra- 
tional Raman  scattering  (RS)  data  for  a hydrogen-air  turbulent  diffusion 
flame  produced  in  a coaxial  jet  combustor.  These  data,  in  the  form  of 
histograms  as  well  as  average  values,  give  radial  profiles  of  tempera- 
ture at  distances  downstream  of  the  fuel  pipe  tip  corresponding  to  1 to 
100  tip  diameters.  Additional  information  on  instantaneous  values  of 
flame  stoichiometry  has  also  been  found.  Experimentation  is  in  progress 
on  coupling  the  RS  apparatus  to  laser  velocimetry  (LV)  equipment,  in 
order  to  produce  essentially  simultaneous  temperature  and  velocity  flow- 
field  information. 


Discussion 

During  the  past  reporting  period,  we  have  focused  upon  preparations 
for  interfacing  the  RS  apparatus  with  LV  instrumentation  for  obtaining 
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turbulent  velocity  data  in  the  same  small  spatial  volume  and  over  the 
same  short  time  duration  as  for  the  RS  data.  We  have  performed  some 
prototype  experiments  to  investigate  key  elements  of  this  work.  We 
have  also  obtained  substantially  more  fluctuation  temperature  values 
from  RS  Stokes/anti-Stokes  data  for  our  turbulent  jet  diffusion  flame 
combustor,  and  are  continuing  adjunct  calibration  and  other  ancillary 
experiments  to  make  the  data  interpretation  more  meaningful. 

The  fluctuation  temperatures  have  been  determined  for  the  co-flowing 
combustor  described  in  previous  reports,  consisting  of  a movable  3 mm 
diameter  fuel  tube  surrounded  by  a 100  mm  diameter  air  pipe.  Turbulent 
diffusion  flames  are  produced  with  a hydrogen  fuel  flow  of  about  585  cm^/s 

3 

and  an  air  flow  of  about  0.1  m /s,  and  the  flame  observed  in  a short  open- 
throat  section.  The  RS  data  have  been  taken  utilizing  a modified  Phase-R 
dye  laser  source,  which  results  in  an  experimental  time  resolution  of 
about  1/2  ps.  The  spatial  resolution,  determined  by  the  collection  optics 

3 

and  the  laser  beam  focusing  properties,  is  roughly  1/2  mm  , with  no  charac- 
teristic dimension  greater  than  about  0.8  mm. 

The  temperature  is  determined  from  the  ratio  R of  integrated  Stokes 
vibrational  scattering  intensity  to  the  integrated  anti-Stokes  vibrational 
scattering  intensity,  viz., 

R = const.  X exp  (hcwo/kT)  (1) 

= const.  X exp  (3374/T)  for  nitrogen. 

where  h is  Planck's  constant,  c is  the  speed  of  light,  wo  is  a vibrational 
term  value  constant  for  the  gas  under  consideration,  k is  Boltzmann's  con- 
stant, and  T is  the  temperature  in  °K.  The  constant  of  proportionality  in 
Equation  (1)  is  best  determined  by  empirical  calibration  of  the  light-scat- 
tering optical  detection  system  through  use  of  a test  medium  of  known 


temperature,  and  has  been  done  here  through  use  of  premixed  hydrogen-air 
flames  burned  on  a porous  plug  burner  for  which  temperatures  were  accu- 
rately determined  in  previous  experiments  (many  of  which  were  carried  out 
under  Project  SQUID  sponsorship). 

The  calibration  experiments  showed  excellent  agreement  among  the 
average  temperatures  of  three  different  flames  (viz.,  from  stoichiometric 
to  lean  premixed  laminar  H2-air  flames),  using  only  one  of  these  flames 
to  determine  the  empirical  constant  in  Equation  (1).  (This  constant, 
which  is  completely  determined  by  the  optical  and  signal  processing  char- 
acteristics of  our  measurement  system  and  is  thus  not  "adjustable"  in  any 
sense,  can  also  be  found  through  use  of  purely  optical  and  electronic  cali- 
bration experiments  for  the  spectroscopic  and  electronic  detection  apparatus. 
Such  measurements  are  now  in  progress.) 

The  relatively  small  statistical  variation  of  these  temperature  data 
was  found  to  be  reasonable  from  the  point  of  view  of  our  experimental  para- 
meters. (For  example,  one  run  of  50  shots  for  which  the  average  value  was 
1807°K  produced  probable  error  bounds  of  ±79°,  corresponding  to  measurement 
bounds  within  which  any  particular  temperature  could  fall  with  50%  proba- 
bility.) Since  the  instantaneous  values  of  temperatures  observed  for  the 
turbulent  coaxial  combustor  were  found  to  have  fluctuations  far  in  excess 
of  the  error  bounds  in  the  laminar  flame  calibration  experiments,  they  are 
convincingly  ascribed  to  the  thermal  fluctuations  of  the  diffusion  flame. 

The  temperature  histograms  (or  probability  distribution  functions)  are 
broader  and  show  more  structure  than  was  expected,  but  the  causes  of  these 
phenomena  are  still  speculative.  The  effect  of  intermittancy  at  the  flame 
boundaries  (i.e.,  in  zones  of  strong  ambient  air  mixing)  is  quite  clear, 
but  additional  histogram-broadening  effects  due  to  flame  oscillations  may 
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exist.  Further  work  in  defining  the  problems  involved  in  the  interpre- 
tation of  the  histograms  is  in  progress;  one  likely  result  will  be  the 
redesign  of  our  combustion  facility  to  incorporate  a closed  throat  square 
test  section  and  a suction  pump  at  the  exhaust  rather  than  a fan  blower. 

Additional  data  of  use  in  interpreting  the  temperature  histograms 
are  given  by  associated  N2  density  histograms,  which  can  be  related  to 
estimates  of  local  flame  stoichiometry.  Current  work  involves  interpre- 
tation of  these  density  histograms  for  their  utility  in  providing  addi- 
tional insight  concerning  the  true  temperature  flowfield. 

The  histograms  obtained  for  our  system  show  strikingly  the  useful- 
ness of  instantaneous  data  acquisition  for  temperature.  Progressing  to 
closely-simul taneous  temperature  and  velocity  data  acquisition  should  re- 
sult in  an  even  clearer  deironstration  of  the  value  of  optical  light- 
scattering diagnostics.  Toward  this  end,  we  are  presently  making  opera- 
tional new  LV  instrumentation  and  are  completing  interface  equipment  for 
coupling  the  RS  and  LV  portions  of  the  experiment.  We  are  also  moving 
from  manual  RS  data  acquisition  (i.e.,  photographic  records  of  oscillo- 
scope traces  for  individually-initiated  experimental  shots)  to  microcom- 
puter-controlled digital  data  acquisition  using  a Nicolet  Explorer  III 
dual  digital  oscilloscope,  which  will  permit  us  to  effectively  couple 
with  the  new  LV  processor.  The  computer  initially  will  use  an  extended 
version  of  the  BASIC  programming  language  to  compute  and  store  tempera- 
ture values  as  well  as  the  velocity  data  from  the  LV  processor.  This 
microcomputer  will  also  be  connected  directly  to  the  GE  CRD  mainframe 
computer  in  order  to  use  effectively  previously-developed  data  reduction 
and  computing  routines. 
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TURBULENCE  MEASUREMENTS 


IN  JETS  FLAMES  AND  COMBUSTORS 


Polytechnic  Institute  of  New  York 
Aerodynamics  Laboratories 


Subcontract  No.  8960-5 
S.  Lederman  - Principal  Investigator 


Introduction 

In  the  last  semiannual  report  of  September  12,  1977,  dat~ 
concerning  concentration  temperature,  velocity,  turbulence 
intensity  and  the  mixedness  parameter  in  a methane  air  CO2 
flame  were  obtained.  In  this  reporting  period,  a particular  set 
of  experiments  concerned  with  the  applicability  of  the  spontan- 
eous Raman  diagnostic  method  to  systems  containing  carbon  mole- 
cules were  conducted.  In  addition  some  experimental  results 
on  a sooty  flame  utilizing  CARS  were  obtained. 

Discussion 

The  applicability  of  the  spontaneous  Raman  effect  to  the 
diagnostics  of  flow  fields  encountered  in  propulsion  and  com- 
bustion has  been  a major  part  in  the  development  of  nonintrusive 

techniques.  As  it  is  well  known,  this  technique  is  applicable 

1-3 

not  only  to  the  measurements  of  concentration  and  temperature  , 

but  also  to  the  measurement  of  fluctuation  intensity  of  the  con- 

4-8 

centration  and  temperature  and  therefore  to  the  fluctuation  in- 
tensity of  the  density  in  a flow  field.  It  has  also  been  shown^, 
that  the  correlation  of  the  second  and  higher  orders  can  also  be 
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measured  using  the  spontaneous  Raman  effect. 

It  is  therefore  important  to  determine  if  the  presence  of 
soot  in  the  flow  field  renders  the  Raman  diagnostic  technique 
useless.  A preliminary  series  of  tests  on  a pure  methane  air  i 

flame,  and  a methane  air  flame  with  the  addition  of  carbon 
particles  were  performed.  Some  results  are  shown  in  Fig.  1 and 

i 

2. 

Using  the  Coherent  Antistokes  Raman  Scattering  system,  1 

some  measurements  were  made  on  a methane  flame  with  varying 
quantities  of  soot.  The  measured  specie  was  the  unburned 
methane  in  the  flame  as  a function  of  X/D.  Some  results  are 
shown  in  Fig.  3 and  4. 

It  should  be  noted  that  an  attempt  was  made  using  the 
Spontaneous  Raman  effect  to  determine  the  unburned  methane  in 
this  flame  using  the  same  Q-switched  Ruby  laser  at  a power  level 
of  over  100  megawatt  peak  power.  No  signal  could  be  detected 
due  to  the  very  low  concentration  of  the  remanent  methane. 
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FJ&-9  CONCENTRATION  OF  UNBURNT  METHANE 


IV.  TURBULENCE 


Semi-Annual  Pronross  Rcoort 


LARGE  SCALE  STRUCTURE  AND  ENTRAINMENT  IN  THE  TURBULENT  MIXING  LAYER 


University  of  Southern  California,  Los  Angeles,  California 
Subcontract  No.  8960-12 


Professor  F.  K.  Browand,  Principal  Investigator 
Mr.  B.  0.  Latigo,  Research  Assistant 


Introduct ion 


Previous  visual  observations  indicate  the  presence  of  large  scale, 
quasi-organized,  vortical  lumps  aligned  across  the  flow  (LSS)  in  the  two 
dimensional  mixing  layer.  The  existence  of  these  structures--docu(nented 
visually  over  a range  of  Reynolds  numbers  extending  from  10^  to  10^  -- 
is  suggestive  of  their  importance  as  a characteristic  feature  of  the 
turbulent  flow.  As  the  mixing  layer  grows  downstream,  the  vortices  must 
necessarily  interact  to  form  larger  vortices.  The  interactions--  to  a 
certain  degree — are  distinct  and  repeatable,  and  it  is  precisely  these 
interactions  which  are  responsible  for  the  growth  of  the  mixing  layer. 
The  presen^  experimental  study,  carried  out  in  a wind  tunnel  at  Reynolds 
numbers  10  , is  intended  to  provide  more  information  about  this  large 
scale  structure. 


D i scussion 


The  procedure  for  studying  large  scale  structure  has  been  to  form  an 
ensemble  average  of  individual  interactions  to  produce  a composite  inter- 
action. The  usefulness  of  this  approach  depends  critically  upon  the 
ability  to  detect  or  sample  repeatably,  similar  interactions.  We  have 
recently  improved  our  sampling  technique  by  introducing  a small  acoustical 
disturbance  at  the  splitter  plate  trailing  edge,  using  a linear  array  of 
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speakers  mounted  in  the  tunnel  roof.  v/ish  to  introduce  a disturbance 
which  will  augment  the  naturally  occurring  interactions  without  altering 
their  fundamental  character. 

The  acoustical  augmentation  produces  several  readily  observable 
effects.  First,  the  large  vortical  structures,  which  grow  dov/nstream 
by  coalescence,  are  phase  locked  with  the  forcing  signal.  The  interactions 
between  these  structures  occur  at  fixed  soatial  locations.  Second,  the 
large  vortical  structures  are  more  organized  and  repeatable  than  without 
forcing.  An  indication  of  this  increased  organization  is  illustrated  in 
figure  1,  where  cross  correlations  between  pairs  of  hot  wires  at  various 
spanwise  separations  are  sho-.n.  The  correlations  on  the  left  result  when 
no  forcing  is  applied.  Now  intermittant  -orcing  is  applied  and  ensembles 
of  velocity  fluctuation  at  different  spanwise  stations  are  formed 
using  the  forcing  to  establish  an  origin  in  time.  Cross  correlations  of 
the  ensemble  averages  are  shown  on  the  right  of  figure  1.  There  is 
practically  no  loss  of  correlation  across  the  span  of  the  wind  tunnel, 
indicating  that  slight  forcing  is  sufficient  to  produce  structures 
which  are  remarkably  two-dimensional. 

There  is  much  work  to  be  done  in  understanding  the  interactions  of 
these  large  scale  structures,  but  preli'^inary  results  suggest  the  following 
concl us  ions ; 

i)  Pairing  interactions  do  occur  at  these  large  Reynolds  numbers, 
and  are  qualitatively  similar  to  the  interactions  at  much  lower  Reynolds 
numbers . 

ii)  The  presence  or  absence  of  the  forcing  does  not  alter  the  major 
features  of  the  interactions. 

iii)  The  interactions  between  the  large  scale  vortices  are  responsible 
for  the  bulk  of  the  observed  Reynolds  stress.  This  result  can  be  seen  in 
two  ways.  First  by  filtering  the  signals,  the  contribution  to  u'v'  from 
different  frequency  bands  can  be  assessed.  Figure  2 shows  the  long-time 
averaged  product  u * v ' for  the  bandwidth  DC  to  2.5  KHz,  compared  to  the 
contribution  from  frequencies  above  280  Hz.  The  passage  frequency  at  this 
downstream  location  is  70  Hz.  It  can  be  concluded  that  the  bulk  of  the 
Reynolds  stress  arises  from  low  frequency  components  comparable  to  the 
passage  frequency.  This  conclusion  is  independent  of  any  assumption  about 
the  existence  of  particular  large  scale  structure.  Some  idea  of  the 
importance  of  the  interaction  of  the  vortical  features  can  be  obtained  by 
calculating  the  contribution  (to  the  long  time  average  u'v')  of  the 
ensemble  averaged  field,  u^Vg  , during  the-  passage  of  the  interaction. 

Figure  2 shows  this  contribution  when  the  interactions  are  assumed  to 
occur  at  the  passage  frequency.  It  is  seen  that  such  interactions  could 
account  for  the  major  portion  of  the  (long  time  averaged)  Reynolds  stress, 
although  the  lower  values  in  the  center  of  the  mixing  layer  may  mean  'we 
arc  not  centered  correctly  in  the  interaction. 

At  present,  two  components  of  velocity  are  being  measured  (with  hot- 
wires) at  a number  closely  soaced  downstream  stations.  At  each  station 
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the  wire  is  traversed  s)ov-/)y  across  the  tlow  to  produce  a continuous  record 
which  is  then  digitized  for  conputet  pi  i.essing. 

We  are  also  preparing  for  publication  some  of  our  earlier  work 
entitled;  "The  Growth  of  the  Tv/o-Dimciis  ional  Turbulent  Mixing  Layer  from 
Turbulent  and  Non-Turbu 1 ent  Initial  Conditions,"  by  F.  K.  Browand  and 
B.  0.  Latigo,  in  preparation. 
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Semi-Annual  Progress  Report 


THE  STRUCTURE  OF  EDDIES  IN  TURBULENT  FLAMES 


University  of  Sheffield,  England 
Subcontract  No . 8960- 30 


Dr.  N.  A.  Chigier,  Principal  Investigator 
Dr.  A.  J.  Yule,  Research  Fellow 


Introduction 


An  experimental  study  is  being  carried  out  on  the  structure 
of  turbulent  propane/air  diffusion  flames.  Within  the  mixing 
regions  of  round  turbulent  jets,  ring  vortices  can  be  detected 
and  followed  as  they  break  down  into  individual  eddies.  The 
flame  characteristics  are  dependent  upon  mixing  of  fuel  and  air 
and  presence  of  coherent  eddies  can  result  in  unmixedness.  The 
main  aim  of  this  study  is  to  relate  high-speed  photographic 
evidence  of  formation  and  movement  of  eddies  and  flamelets  with 
time-resolved  measurements  of  velocity  and  temperature. 

Variation  of  velocity  with  time  is  measured  using  a single 
particle  counting  forwardscatter  laser  anemometer,  specially 
developed  for  conditional  sampling  of  velocity  in  turbulent  flames. 
Temperature  variation  with  time  is  measured  with  fine  wire  thermo- 
couples with  direct  con^sensation  by  interfacing  with  a computer. 


Discussion 


The  apparatus,  consisting  of  a central  burner  installed  in  a 
settling  chamber,  has  been  constructed  to  provide  uniform  constant 
velocity  with  varying  turbulence  levels  in  the  central  jet  and  in 
the  surrounding  secondary  flow.  Propane  is  partially  premixed 
with  air  containing  Si02  seeding  particles  for  supply  to  the 
central  burner.  A cyclone-seeding  apparatus  has  been  constructed 
for  the  central  flow  seeding  and  the  secondary  flow  is  seeded 
separately  by  a high  throughput  fluidized  bed  seeder.  In  order  to 
minimize  biasing  errors  in  velocity  measurements  in  the  flame. 
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seeding  levels  in  the  inner  and  outer  flows  require  to  be  equal 
and  maintained  at  a constant  level  throughout  the  experiment. 
Problems  have  been  encountered  with  seeding  particles  settling  on 
screens  introduced  for  smoothing  the  flow.  The  laser  anemometer 
has  been  used  to  measure  exit  velocity  profiles  and  achieve 
optimum  design  of  the  screens  in  order  to  produce  low  turbulence 
profiles  at  the  exit  of  both  the  inner  burner  flow  and  the  outer 
air  flow. 

Photographic  studies  have  been  made  of  partially  premixed 
propane/air  diffusion  flames.  These  studies  have  been  used  to 
locate  the  flame  region  and  select  flow  parameters  for  detailed 
measurement  studies.  As  in  the  case  of  cold  vortex  studies 
previously  made  by  Yule  at  Southampton,  ring  vortices  can  be 
clearly  identified  in  the  flames  for  a wide  range  of  Reynolds 
number  up  to  10^.  It  was  found  that  the  existence  of  a vortex 
ring  region  depended  critically  on  the  initial  turbulence  level  at 
the  exit  from  the  burner . When  this  turbulence  intensity 
increased  above  6% , the  vortex  rings  were  observed  to  lose  their 
symnetry  very  rapidly. 

During  experimental  runs,  it  was  found  that  the  high  seeding 
densities  necessary  for  high  frequency  response  velocity  measure- 
ments resulted  in  rapid  clogging  of  fine  gauzes  installed  for  flow 
smoothing.  In  order  that  the  initial  conditions  could  remain 
constant,  with  low  turbulence  levels,  it  was  necessary  to  frequently 
dismantle  and  clean  the  burner  assembly.  It  appears  that  there  is 
an  incompatibility  in  achieving  low  turbulence  levels  using  fine 
wire  gauzes  and  the  use  of  particles  for  seeding.  A redesigned 
burner  is  now  under  construction,  which  is  intended  to  overcome 
this  difficulty. 

During  the  next  six  months,  vortex  ring  flames  will  be  studied 
by  simultaneous  measiarement  of  velocity  using  the  laser  anemometer, 
and  temperature  using  the  fine  wire  thermocouples.  The  initial 
series  of  measurements  will  be  made  in  the  flame  to  obtain  mean, 
rms,  probability  density  function  and  cross-correlations  with  time 
delay  for  velocity  and  temperature.  Visual  examination  of  time 
histories  and  photographic  studies  - with  temperature  simultaneously 
recorded  on  the  film  - will  be  used  as  a basis  for  designing  a 
conditional  sampling  technique,  which  will  then  be  used  to  study 
the  structure  of  individual  vortex  rings. 
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INTRODUCTION 


It  is  the  purpose  of  our  research  to  illuminate,  and  to 
provide  quantitative  measures  of  the  fundamental  phenomena 
which  are  responsible  for  the  strongly  enhanced  molecular 
diffusivity  effects  in  a turbulent  mixing  field.  The  presence 
of  these  effects  is  of  obvious  importance  in  the  combustion 
process;  their  full  exploitation  requires  an  understanding  of 
their  dependence  upon  the  character  of  the  turbulence  field. 

One  approach  toward  this  understanding  is  to  examine  the 
results  of  controlled  variations  in  the  governing  parameters 
of  experiments  which  are  (1)  sufficiently  simple  that  the 
cause/effect  relationships  are  least  ambiguous  and  (2) 
sufficiently  similar  to  the  technological  problem  that  the 
phenomena  of  interest  are  preserved.  Our  experiments  examine 
the  mixing  of  two  distinct  rectangular  volumes  by  light  scattering 
measurements  from  the  central  region  of  a closed  mixing  chamber. 
The  nature  of  the  experimental  facility  allows  the  initial 
turbulence  structure  in  the  two  volumes  to  be  individually 
controlled  and  stable,  unstable  or  neutral  density  mixing  may 
be  investigated. 

The  first  set  of  experiments  using  this  unique  facility 
has  been  completed.  The  lower  volume  was  contaminated  with 
submicron  particles,  which  allowed  the  non-dif f usive  mixing  of 
the  upper  and  lower  volumes  to  be  observed  by  noting  the 
presence  or  absence  of  the  marking  contaminant  in  small  sample 
volumes.  The  volumes  are  defined  by  a 0.25  mm  (focused)  laser 
beam  and  the  0.35  mm  height  was  defined  by  a helical  slit  in  a 
rotating  disc  located  at  the  plane  of  tne  recreated  laser  beam. 
Each  scan  (by  a slit)  sampled  205  locations  along  a 74  mm 
line  parallel  to  the  concentration  gradient;  the  scan  was  executed 
in  3.11  msec  (=  0.768  M/U^)  and  the  concentration  r(Z,t)  was 
denoted  as  1 or  0 depending  upon  the  presence  or  absence  of  the 
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contaminant.  A data  base  of  69  experiments  with  air/air  and 
100  experiments  with  Freon  12/Freon  12  was  acquired.  These 
experiments  were  to  establish  the  experimental  technique  for 
the  simplest  case  of  neutral  density  mixing  and  they  were  to 
establish  the  influence  of  the  kinematic  viscosity  on  the  mixing 
field  = 5.8  12^’  These  objectives  have  been  met. 

The  experiments  and  the  associated  analyses  form  the  basis  of 
the  Ph.D.  thesis  for  Mr.  K.  C.  Cornelius.  The  thesis  was 
defended  in  February  and  is  in  the  stage  of  final  revisions  and 
typing.  The  SQUID  report  to  be  written  from  it  will  be  forth- 
coming . 

DISCUSSION 


The  objective  reporting  of  the  experiments,  which  was  noted 
as  being  essentially  complete  in  the  previous  semi-annual 
progress  report,  has  been  considerably  extended  by  the  develop- 
ment of  additional  theoretical  analyses.  Width  and  displacement 
measures  have  been  obtained  for  the  ensemble  average  of  the 
instantaneous  scans  once  these  are  shifted  into  alignment  with 
respect  to  their  instantaneous  centers.  The  measures,  so 
defined,  reveal  that  the  centerline  displacements  are  insensitive 
to  the  six  fold  change  in  kinematic  viscosity  (as  would  be  expec- 
ted but  that  the  relative  width  of  the  instantaneous  scan  is  de- 
pendent upon  this  fluid  property.  (Note,  U and  M are  identical 
for  the  air  and  Freon  12  test  gases;  v . =°5.8  12^' 

is  also  shown  that  the  ensemble  average^^iixing  figf§°aii€h  is 
closely  approximated  by  the  same  measure  from  the  instantaneous 
scans  and  that  the  Reynolds  number  independent  gross  transport 
makes  only  a small  contribution  to  the  observed  width. 

The  second  theoretical  development  exploits  the  observation 
that  — <r(Z,t)>  (the  ensemble  averaged  concentration  distribution) 
can  be  directly  related  to  the  dispersion  of  a single  particle  — 
in  order  to  infer  the  Lagrangian  microscale  and  macroscale. 

The  presence  of  the  decaying  turbulence  field  is  acknowledged 
by  utilizing  the  basic  proposal  of  Townsend  (1954)  and  extending 
the  analysis  of  Shlien  and  Corrsin  (1974).  This  improved  analysis 
allows  the  data  to  be  evaluated  in  terms  of  the  laboratory 
(normalized)  time  coordinate  instead  of  in  the  stretched  time 
coordinate;  the  Shlien  and  Corrsin  data  are  reinterpreted  with 
the  new  formulation.  With  respect  to  the  present  results,  the 
new  theory  supports  the  observation  that  the  air  dispersion 
is  greater  than  that  for  the  Freon.  The  difference  exists  in 
the  Reynolds  number  dependence  of  the  Lagrangian  micro  scale; 
details  are  provided  in  the  full  report. 

The  identification  of  the  kinematic  viscosity  influence 
on  the  mixing  (dispersion)  phenomenon  is  important  for  the 
projected  large  density  difference  experiments.  It  was  earlier 
suspected  that  the  differences  in  the  Freon  12  and  air  results 
were  a result  of  sub-scattering  volume  mixing.  The  results 
from  the  instantaneous  scans  and  the  dispersion  measurements 
provide  independent  confirmation  that  this  potential  problem 
was  not  the  basis  for  the  observed  differences.  The  projected 
experiments  will  have  to  be  interpreted  with  the  understanding 
that  the  Lagrangian  microscale  is  larger  in  the  more  viscous 
gas . 
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Introduction 

This  research  addresses  problems  related  to  the  turbulent 
heterogeneous  flows  which  arise  in  a variety  of  propulsive  devices 
when  reactants  and  products  mix  and  react.  The  effort  is  both  experi- 
mental and  theoretical;  the  experimental  program  concerns  exploitation 
and  extension  of  the  multiple  sensor  "hotwire"  technique  of  Way  and 
Libby  which  permits  time-resolved  and  space-resolved  measurements 
of  velocity  and  concentration  of  one  light  species,  e.  g. , helium,  in  a 
mixture  of  light  and  heavy  gases  under  isothermal  conditions.  The 
application  of  this  technique  in  the  present  research  is  to  a confined 
internal  flow  corresponding  to  an  idealized  combustor.  The  related 
theoretical  work  supports  the  experimental  effort  and  attempts  to  extend 
the  results  thereof  to  flow  situations  of  more  practical  concern,  e.g.,  to 
chemically  reacting  flows. 

Discussion 

In  our  previous  Semi-Annual  Progress  Report  we  described  work 
underway  on  premixed  combustion.  We  have  for  several  years  been 
collaborating  with  Professor  K.  N.  C.  Bray  and  Dr.  J.  B.  Moss  of  the 
University  of  Southampton  on  applications  of  the  Bray-Moss  model  of 
premixed  combustion  to  planar  turbulent  flames,  both  oblique  and  normal. 
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Our  most  recent  effort  has  been  concerned  primarily  with  the  influence 
of  large  but  finite  values  for  the  two  non-dimensional  parameters 
determining  the  behavior  of  such  flames,  a turbulence  Reynolds  number 
and  a Damkohler  number.  The  original  theory,  termed  first-order, 
assumes  that  these  two  parameters  are  infinite  and  we  wished  to  examine 
the  effect  of  slight  deviations  from  this  limit  in  order  to  expand  the  range 
of  comparison  with  data.  This  work  has  been  completed  and  submitted 
for  possible  publication  as  indicated  by  reference  1. 

In  the  course  of  the  study  leading  to  reference  1 considerable 
attention  was  devoted  to  comparing  the  predictions  of  the  first-order 
calculations  to  experimental  data  for  normal  flames.  In  previous  work 
we  had  focused  on  such  comparison  for  highly  oblique  flames  and  had 
established  that  a value  of  the  one  empirical  parameter  in  the  first-order 
theory,  namely  <I>,  equal  to  0.  I was  required  to  obtain  agreement.  The 
experimental  data  for  normal  flames  at  the  highest  turbulence  Reynolds 
numbers  involve  considerable  scatter  so  that  an  unequivocal  value  for 
4*  is  difficult  to  establish.  Nevertheless,  a value  of  1.4,  an  order  of 
magnitude  greater  than  that  needed  to  achieve  agreement  for  highly 
oblique  flames,  seems  indicated.  We  thus  observe  another  example  of 
the  shattered  hope  for  "universality"  in  the  constants  appearing  in  turbulence 
phenomenology.  There  are,  of  course,  a variety  of  possible  explanations 
for  the  discrepancy;  in  terms  of  the  theory  the  use  of  gradient  assumptions 
under  the  conditions  of  high  rates  of  strain  such  as  arise  in  turbulent 
flames  and  the  neglect  of  the  effect  of  weak  pressure  gradients  when 
large  density  variations  are  present  are  two  candidates  for  shortcomings 
in  the  theory.  With  respect  to  experiment  the  wide  scatter  in  available 
data  suggests  that  the  experimentalists  cannot  agree  on  a fundamental 
quantity  of  concern  to  the  theoretician,  the  turbulent  flame  speed.  We 
continue  research  on  this  problem. 

During  the  past  few  months  a new  avenue  of  theoretical  work  related 
to  premixed  combustion  has  been  initiated.  We  have  been  examining  the 
influence  of  acceleration  as  may  arise  in  unsteady  turbulent  flames  which 
are  driven  hydrodynamically  or  in  horizontal  flames  propagating  steadily 
in  either  the  upward  or  downward  direction.  There  appears  to  have  been 
no  previous  examination  of  this  influence.  Suggestions  are  made  in  the 
literature  of  Taylor  instability  as  a cause  for  the  onset  of  turbulence  in 
accelerated  laminar  flames  but  we  are  concerned  with  the  alteration  of 
turbulent  flame  behavior  as  a result  of  acceleration. 

From  an  examination  of  the  describing  equations  and  from  a 
consideration  of  the  only  related  literature,  that  concerned  with  turbulence 
involving  buoyancy,  we  have  deduced  that  acceleration  can  have  two 
effects;  the  first  is  the  alteration  of  the  turbulent  kinetic  energy  balance 
by  a source  term  proportional  to  acceleration  and  to  the  extent  of  heat 
release.  The  second  effect  is  a more  subtle  one,  namely  the  alteration 
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of  the  flux  of  product  which  in  Favre-averaging  is  pu”c*  . In  our 
previous  calculations  related  to  unaccelerated  flames  this  flux  is  modeled 
in  terms  of  a gradient  approximation.  Inclusion  of  the  effect  of  acceleration 
on  this  flux  appears  to  require  a new  theory  which  would  apply  to  both 
steady  and  accelerated  flames.  Such  a theory  for  normal  flames  has  been 
formulated  but  no  results  are  at  present  available. 

Our  first  calculations  of  the  effect  of  acceleration  on  normal  flames 
were  based  on  alteration  of  only  the  turbulent  kinetic  energy  balance  and 
show  that  acceleration  can  lead  to  a significant  increase  in  the  turbulent 
intensity  through  the  flame.  This  is  in  contrast  to  the  reduction  in  that 
energy  due  to  dilatation  when  the  flame  is  unaccelerated.  There  is  also 
an  increase  in  turbulent  flame  speed.  However,  in  view  of  the  possible 
importance  of  acceleration  on  the  flux  of  product  we  do  not  consider  these 
results  definitive  and  await  the  completion  and  assessment  of  the  new 
theory  before  making  a definite  conclusion  on  the  behavior  of  flames 
undergoing  acceleration. 

During  the  period  covered  by  this  Report  we  have  completed  for 
a DOE  Workshop  a critique  of  problems  in  the  phenomenology  of  turbulent 
reacting  flows  (reference  2)  and  have  submitted  it  for  publication  in  the 
proceedings  of  that  Workshop. 

Preparations  for  our  next  experiment  involving  helium-air  mixing 
are  nearing  fruition  and  tests  will  be  run  in  the  next  several  weeks. 

Involved  will  be  two  distinct  experiments.  The  first  will  correspond 
to  the  tests  reported  in  reference  3 and  will  involve  helium  at  ambient 
temperature  discharged  into  a moving  airstream.  These  tests  will 
complement  and  supplement  the  data  in  reference  3 and  will  permit 
that  research  to  be  submitted  for  publication  in  the  open  literature.  The 
second  experiment  will  involve  heating  of  the  helium  in  order  to  provide 
data  on  the  behavior  of  two  scalars  in  turbulent  shear  flows.  There  are 
at  the  present  no  data  of  this  sort  and  yet  there  is  considerable  fundamental 
interest  since  the  molecular  properties  of  the  two  scalars  differ  significantly. 
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Introduction 


The  objective  of  this  research  is  to  obtain  a better  understanding 
of  the  turbulent  mixing  processes  that  occur  in  mixing  layers  between 
gas  streams  of  different  velocities  and  densities.  Such  mixing  layers 
are  often  a basic  element  in  flows  which  occur  in  propulsive  devices; 
examples  of  problems  to  which  the  research  is  relevant  include  tur- 
bulent combustion,  jet  noise,  and  thrust  augmentation.  The  research 
has  proctreded  along  two  parallel  lines.  On  the  one  hand,  we  have 
been  making  measurements  of  various  statistical  properties  of  the 
mixing  region  and  their  dependence  on  parameters  such  as  Reynolds 
number,  velocity  ratio  and  density  ratio.  Such  information  provides 
important  inputs  for  engineering  models  and  calculation  methods. 

On  the  other  hand,  we  have  been  using  the  quantitative  measurements, 
e.  g.  , time-  and  space  - resolved  concentration  measurements,  together 
with  flow  visualization  to  identify  and  describe  the  physical  processes 
occuring  in  such  mixing  regions.  Better  understanding  of  the  physics 
is  important  for  the  development  of  more  realistic  computing  models 
and  also  for  suggesting  how  turbulent  mixing  might  be  controlled 
or  modified. 


Discussion 

The  main  effort  continues  to  be  devoted  to  the  question  of  how 
three  dimensionality  develops  in  plane  mixing  layers  which,  from 
previous  work  in  this  and  other  laboratories,  are  now  known  to  be 
dominated  by  large  vortices  or  rollers  which  have  fairly  good  span- 
wise  organization  (Refs.  1,  2).  In  recent  work  under  this  contract 
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(Ref.  3),  the  presence  of  streamwise  streaks  was  discovered  on 
spark  shadow  or  schlieren  views  taken  normal  to  the  plane  of  the 
mixing  layer.  An  example  of  a schlieren  picture  is  shown  here 
in  Figure  1 which  shows  simultaneous  side  and  plane  views  of  a 
mixing  layer  between  streams  of  helium  and  nitrogen  at  velocities 
752  and  227  cm/sec.,  respectively,  at  a pressure  of  4 atm.  In 
Ref.  3 it  was  suggested  that  these  streaks  mark  the  edges  of 
streamwise  vortices  which  are  created  by  some  kind  of  seconday 
instability  on  the  main  vortices,  possibly  of  the  Gortler  or  Taylor 
type.  We  are  engaged  in  a series  of  experiments  to  further  elu- 
cidate the  nature  of  these  secondary  organized  structures  and  to 
relate  them  to  various  parameters  of  the  mixing  layer.  We  believe 
that  their  appearance  and  subsequent  development  is  an  important 
link  in  the  development  of  the  three  dimensionality  and  small  scales 
that  must  be  present  in  turbulent  flow  at  high  Reynolds  number. 

They  are  also  thought  to  be  related,  in  some  way,  to  the  enhanced 
internal  mixing  which  occurs  above  a certain  Reynolds  number, 
as  found  in  Ref.  3. 

On  repeated  pictures  like  that  in  Figure  1,  it  was  noted  that 
the  streaks  appeared  to  occur  in  more  or  less  the  same  spanwise 
positions.  This  was  confirmed  by  making  long  time  exposures, 
an  example  of  which  is  shown  in  Figure  2a  for  the  same  conditions 
as  in  Figure  1.  Most  of  the  recent  work  has  been  centered  on 
studying  these  long-time  exposure  patterns  and  their  changes  with 
various  parameters. 

One  question  is  whether  the  streaks  originate  from  fixed 
disturbance  points  upstream  in  the  settling  chamber  or  in  the 
boundary  layers  of  the  splitter  plate,  or  whether  their  spacing  is 
govi-rned  by  an  instability  criterion.  Accordingly,  various  changes 
weri!  made  in  the  screen  configuration  in  the  settling  chamber,  in 
the  length  of  the  splitter  plate  and  in  the  trailing  edge  of  the  latter. 
Also  disturbances  in  the  nature  of  small  vortex  generators  were 
deliberately  introduced  in  the  entrance  flow.  Secondly,  the  various 
parameters  of  the  flow,  in  particular  the  Reynolds  number,  were 
varied.  From  all  these  we  tentatively  conclude  that  the  mean 
pattern  of  the  streaks  is  not  tied  to  upstream  disturbances  but  is 
determined  by  flow  criteria.  Some  of  the  results  on  which  we  base 
this  conclusion  are  the  following. 

1.  The  pattern  of  streaks  depends  on  Reynolds  number  Re. 
When  velocity  and  pressure  are  changed  in  such  a way  as  to  keep 
Re  constant,  the  pattern  remains  unchanged.  When  Re  is  increased 
the  pattern  shifts  upstream  and  the  average  spacing  decreases,  and 
conversely  for  decrease  in  Re.  If  we  denote  the  position  of  the 
front  of  the  pattern  by  x^  then  the  dependence  on  Reynolds  number 
Re  = (U^-U2)(x^-Xq)  / for  various  values  of  the  velocity  ratio 

U2/U1  is  as  shown  in  the  following  table,  where  x^  is  the  iiffective 
origin  of  the  shear  layer  and  is  the  kinematic  viscosity  on  the 

(lower  speed)  nitrogen  side. 
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0.  14 


0.  30 


0.  38 


0.  50 


»2/^l 

10''^R<r  5.9  3.9  -5.1  1.5 

2.  The  i-ncan  spacing  of  the  streaks  or  banris  scalt^s  with 
the  momentum  thickness  of  the  initial,  laminar  shear  layer  approxi- 
mately as  X = 18F).  The  adjustments  in  tht'  spacing,  when  flow  con- 
ditions are  changed,  is  by  appearance  or  disappearance  of  some  of 
the  bands  and  slight  displacement  of  some  of  the:  others. 

3.  The  spanwise  positions,  but  not  the  mean  spacing  of 
the  bands,  could  be  altered  by  changing  the  screen  configuration  in 
the  settling  chamber. 

Still  another  question  concerns  the  relation  of  these  patterns 
to  the  enhanced  mixing  which  was  noted  in  Ref,  3.  In  figure  2b 
the  streamwise  development  of  a mixedness  parameter,  which  was 
measured  in  Ref.  3,  is  scaled  to  the  flow  picture  in  i'igure  2a. 

From  this  it  appears  that  the  increase  in  mixedness  occurs  toward 
the  downstream  side  of  the  streak  pattern.  Assuming  that  these 
streaks  correspond  to  streamwise  vortices,  it  may  be  conjectured 
that  the  enhanced  mixing  corresponds  to  some  further  instability 
or  small  scale  breakdown  of  those  vortices.  In  Ref.  3 it  was  also 
found  that  the  enhanced  mixedness  was  accompanied  by  an  increase 
in  power  spectral  content  at  high  values  of  wave  number.  Those 
wave  numbers  are  considerably  higher  than  what  would  be  inferred 
from  the  mean  streak  spacing,  again  suggesting  a breakdown  to 
smaller  scales.  New  experiments  designed  to  correlate  various 
details  of  the  fluctuating  concentration  field  and  the  streak  patterns 
are  underway. 


Notes  and  References 


1.  Brown,  G.  L.  and  Roshko,  A.  1974  On  density  effects  and 
'large  structure  in  turbulent  mixing  layers.  J.  of  Fluid  Mechanics 

775-816. 

2.  Winant,  C.  D.  Browand,  F.  K.  1974  Vortex  pairing;  the 
mechanism  of  turbulent  mixing-layer  growth  at  moderate  Reynolds 
number.  J.  of  Fluid  Mechanics  63,  237-255. 

3.  Konrad,  John  Harrison  1976  An  experimental  investigation  of 
mixing  in  two-dimensional  turbulent  shear  flows  with  applications 
to  diffusion-limited  chemical  reactions.  Project  SQUID  Technical 
Report  CIT-8-PU. 
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Professor  M.  S.  Uberoi,  Principal  Investigator 


Discussion 


Role  of  axial  flow  in  laminar  and  turbulent  swirling  flows.  Tn  an 
extensive  theoretical  investigation  (Project  SQUID  Technical  Report 
UC-l-PU,  Mechanisms  of  Decay  of  Laminar  and  Turbulent  Vortices,  1977) 
we  have  shown  that  the  variation  of  axial  flow  plays  a very  Important 
role  in  the  development  of  swirling  flows.  We  are  investigating  methods 
of  calculating  these  axial  flows  in  some  swirling  flows  of  practical 
Interest.  Our  attention  is  first  confined  to  the  case  of  a potential 
vortex  with  a concentrated  distribution  of  vortlcity  on  the  axis  of  the 
vortex. 


Coherent  structure  and  Instabilities  in  the  initial  region  of  a plane 
turbulent  air  jet  of  variable  initial  temperature.  Instabilities  and 
coherent  structure  of  isothermal  and  heated  air  jets  are  examined  using 
shadowgraph,  velocity  and  temperature  probes  in  Reynolds  number  range 
10  - 10  based  on  jet  width.  In  the  two  mixing  regions  at  the  jet  edges 

shadowgraphs  show  a very  coherent  structure  of  two  symmetric  spatially 
growing  eddies,  which  engulf  the  whole  jet  and  then  break  up  to  produce  a 
turbulent  jet.  The  distance  required  for  the  roll  up  of  the  vortex  sheet 
or  the  development  of  eddies  and  their  breakdown  depends  on  the  initial 
velocity  profile  of  the  jet.  The  spectra  of  velocity  and  temperature  fluc- 
tuations at  a point  do  not  show  a distinct  frequency,  there  being  considerable 
jitter,  although  a central  frequency  can  be  assigned  to  the  spectra. 

Strouhal  number  based  on  the  central  frequency  was  measured  as  a function 
of  jet  Reynolds  number.  The  roll  up  of  the  two  vortex  sheets  at  the  jet 
edges  into  distinct  eddies  decreased  as  the  Influence  gf  the  buoyancy  was 
Increased  by  Increasing  initial  jet  temperature  to  300  C. 


We  are 
preparing  a 


in  the  process  of  completing  one  phase  of  this  work  and  are 
report  for  publication. 
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SECOND-ORDER  CLOSURE  MODELING 
OF  TURBULENT  COMBUSTION 


Aeronautical  Research  Associates  of  Princeton,  Inc. 
Princeton,  New  Jersey 
Subcontract  No.  8960-26 


Ashok  K.  Varma,  Principal  Investigator 
Coleman  duP.  Donaldson 


Introduction 


The  interaction  between  turbulence  and  chemistry  is  of  considerable 
importance  in  determining  combustion  efficiency  and  pollutant  formation 
as  well  as  other  combustion  characteristics  in  many  combustion  and  pro- 
pulsion systems.  This  research  prc)gram  is  directed  towards  the  study  of 
the  turbulence-chemistry  interaction  using  a complete  second-order 
closure  approach. 

A second-order  closure  procedure  re<}viires  the  development  of  closure 
models  for  the  higher-order  correlations  that  appear  in  the  transport 
equations  for  the  means  and  the  second -C)rder  correlations.  We  have  pro- 
posed a delta  function  "typical  eddy"  model  for  the  joint  probability 
density  function  for  all  the  scalars.  We  believe  that  a new  principle 
has  been  formulated  which  permits  rational  closures  for  turbulent  mixing 
and  turbulent  reacting  flows. 


Discussion 


Recently  a major  advance  has  been  made  In  the  procedure  for  the 
construction  of  delta  function  typical  eddy  models.  In  the  past  we 
have  struggled  with  defining  a criterion  for  selecting  the  proportions 
of  the  species  in  the  mixed  cells  ( In  a two  species  model  and 
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, 05  and  6^  in  a three  species  model).  In  the  previous  progress 

report  (Ref.  1)  we  had  proposed  using  the  midrange  value  of  the  s^ 
moment  to  select  03.  However,  we  now  believe  that  the  following  cri- 
terion should  be  used  for  determining  the  typical  eddy  structure. 

Minimum  Entropy  of  Mixing  Principle.  The  free  parameters  in  the 
delj; a function  typic al_  eddy  mod e 1 should  be  chosen  to  minimize  the 
^nt ropy  o f mixing . 

The  transport  equations  for  the  means  and  various  second-order 
correlations  Incorporate  the  entropy-increasing  effects  of  molecular 
transport  processes,  and  this  is  reflected  in  the  values  of  the  means 
and  correlations  at  any  given  point  in  the  flow.  It  is  important  that 
the  modeling  procedure  not  add  any  more  entropy  to  the  system,  that  is, 
there  should  be  no  entropy  of  modeling  or  entropy  increase  due  to  our 
ignorance.  This  can  be  accomplished  by  selecting  the  model  structure, 
within  the  allowable  range  of  parameters,  that  has  the  least  entropy. 

The  notion  can  also  be  understood  from  the  viewpoint  of  large  eddy 
structures.  Consider  the  mixing  of  two  initially  separate  streams  of 
inert  gases  a and  3 • The  large  eddies  in  a turbulent  flow  cause 
stretching  and  folding  and  Intertwining  of  the  streams,  but  do  not  cause 
an  Increase  in  the  entropy  of  the  system.  The  entropy  of  mixing  reflects 
only  the  molecular  mixing  processes  and  therefore,  selecting  the  minimum 
entropy  model  is  consistent  with  the  large  structures  viewpoint  of  turbu- 
lent flow. 

For  a two  species  model,  the  parameter  03  is  selected  on  the  basis 
of  minimum  entropy.  For  a three  species  model,  one  has  to  simultaneously 
determine  , 05  and  to  minimize  the  entropy.  In  practice,  this 

concept  is  very  convenient  to  use  and  is  working  well.  We  have  carried 
out  tests  on  a two  species  model  in  the  past  month  and  some  typical  results 
are  discussed  below.  In  practice,  the  entropy  appears  to  be  a smoothly- 
varying  function  of  the  parameter  03  and  the  selection  of  03  for 
minimum  entropy  is  quite  straightforward.  A few  examples  of  the  minimum 
entropy  "typical  eddy"  structures  at  various  points  in  the  statistically 
valid  moment  space  are  shown  below. 


The  model  structure  is  computed  as  a function  of  03  for  specified 
values  of  the  moments  ci  , ^ , so^  and  s^^  . In  general,  valid  model 

structures  ^te  obtained  for  a range  of  values  of 

03  . The  program  calculates  the  total  entropy  of  the  structures  and  se- 
lects the  value  of  03  for  minimum  entropy. 
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Model  Structure 
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Consider  tlie  point  in  moment  space. 
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Further  tests  are  continuing  to  ascertain  whetlier  the  minimum 
entropy  structure  is  always  a three  cell  structure  as  was  the  case 
in  the  above  two  examples,  and  if  this  were  so,  to  attempt  to  formulate 
even  simpler  rules  to  determine  the  structure. 

We  have  also  made  progress  in  other  areas  of  our  effort. 

Second-Order  _C2osi^e_  Program  Development 

The  second-order  closure  program  including  the  density  correlation 
equations  is  now  operational.  The  equations  for  p'a'  , o'S'  and  p’y' 
are  solved  in  pass  11  and  the  equations  for  p *h*  and  p ' p ’ are  solved 
in  sequence  in  pass  12  of  the  RSL  computer  program. 
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The  only  remaining  task  now  is  to  incorporate  the  solution  procedure 
for  the  two  species  typical  eddy  model  as  a subroutine  to  the  RSL  program. 
On  the  completion  of  this  phase  of  the  effort,  we  will  compare  the  results 
of  second-order  closure  calculations  for  a two  species  mixing  flow  to 
experimental  measurements  of  Brown  and  Roshko  and  Konrad. 

Work  is  also  continuing  on  the  development  of  the  correct  statistical 
constraints  for  a three  species  system.  Tests  can  then  be  carried  out  to 
apply  the  minimum  entropy  principle  to  a three  species  "typical  eddy" 
model . 
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APPENDIX  A 
TECHNICAL  REPORTS 
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